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Abstract 

We present a definition of the four-dimensional helicity (FDH) regularization scheme valid for 
two or more loops. This scheme was previously defined and utilized at one loop. It amounts to a 
variation on the standard 't Hooft-Vcltman scheme and is designed to be compatible with the use 
of helicity states for "observed" particles. It is similar to dimensional reduction in that it main- 
tains an equal number of bosonic and fermionic states, as required for preserving supersymmctry. 
Supersymmetry Ward identities relate different helicity amplitudes in supersymmetric theories. 
As a check that the FDH scheme preserves supersymmetry, at least through two loops, we explic- 
itly verify a number of these identities for gluon-gluon scattering (gg — > gg) in supersymmetric 
QCD. These results also cross-check recent non-trivial two-loop calculations in ordinary QCD. 
Finally, we compute the two-loop shift between the FDH coupling and the standard MS coupling, 
a s . The FDH shift is identical to the one for dimensional reduction. The two-loop coupling shifts 
are then used to obtain the three-loop QCD j3 function in the FDH and dimensional reduction 
schemes. 
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1 Introduction 



It is of great utility to regularize the divergences of quantum field theory in a way that manifestly 
preserves the symmetries of the theory. The most widely utilized technique for preserving gauge sym- 
metry is dimensional regularization, which can simultaneously handle the ultraviolet and infrared 
divergences of massless gauge theory. There are actually several common variants of dimensional 
regularization, differing in their treatment of the Lorentz vector indices associated with gauge parti- 
cles. The different variants have advantages and disadvantages, depending on the application. The 
conventional dimensional regularization scheme (CDR) [Qj is conceptually the simplest and most 
widely used variant. In this scheme one uniformly continues all momenta and vector polarizations to 
D = 4 — 2e dimensions. The 't Hooft-Veltman (HV) scheme [Q] is closely related, differing only in the 
treatment of "observed" states, which remain in four dimensions. The HV scheme is especially con- 
venient for computing helicity amplitudes — as has proven very useful in QCD calculations (3|, [|, ||, ||] 
- because only the amplitudes with the four-dimensional helicity values ±1 (or 0, for massive vector 
bosons) have to be evaluated. For supersymmetric theories, the most commonly used regularization 
scheme is Siegel's dimensional reduction (DR) scheme J/]], since it preserves super symmetry. This 
scheme has also found use in certain non-super symmetric calculations ||, |J. Other useful schemes 
for supersymmetric theories include the holomorphic and NSVZ schemes ]l0| . 

Another scheme that has been used at one loop is the four-dimensional helicity (FDH) scheme ||. 
The FDH scheme seeks to combine the natural use of helicity states (as in the HV scheme) with 
the preservation of supersymmetry. It is similar to the DR scheme in maintaining the number of 
physical states at their four-dimensional values. In the FDH scheme, however, the algebraic rules 
are defined with the notion that D > 4, but with an analytic continuation to bring the number of 
physical states back to their D = 4 values. This may be contrasted with the DR scheme, where one 
takes D < 4, viewed as a dimensional compactification. The distinction between D < 4 and D > 4 
is relevant only when an explicit basis of external states is required for a calculation; inside the loops 
in either case all indices on fields are treated as four-dimensional. At one loop, the relationship 
between the FDH and DR schemes has been previously discussed in ref. 

Preserving supersymmetry is useful even for higher order calculations in QCD. Although QCD 
itself is not a supersymmetric theory, one can slightly modify QCD by altering the color repre- 
sentations and multiplicities so that it becomes supersymmetric. The amplitudes in such a theory 
are closely related to those in QCD, yet they must satisfy non-trivial supersymmetry Ward identi- 
ties |l2|, 13], as long as the regulator preserves supersymmetry. These identities can therefore provide 
an independent means for checking a non-trivial QCD calculation. The supersymmetry identities 
on scattering amplitudes are phrased in terms of the helicity basis; thus the FDH scheme is ideal 
for using supersymmetry in this way. Of course one would also want to have available a simple 
way to convert amplitudes computed in the standard variants of dimensional regularization, such as 
the 't Hooft-Veltman (HV) or conventional dimensional regularization (CDR) schemes, into those 
computed in the FDH scheme, or vice-versa. 

Though we do not have a proof that the FDH scheme preserves supersymmetry to all orders of 
perturbation theory, it inherently maintains the number of physical states at their four-dimensional 
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values. Moreover, from explicit calculations the scheme is known to preserve supersymmetry at one 



loop [11]. In this paper we shall explicitly verify its preservation at two loops, by checking various 
supersymmetry identities involving gg — * gg helicity amplitudes. 

The issue of higher-loop regularization is timely in light of the recent substantial progress in 
the calculation of two-loop scattering amplitudes. Until recently, no such amplitudes depending on 
more than a single kinematic variable were known. Now several such computations have appeared. 
The first of these were gluon-gluon scattering amplitudes in the special cases of maximal (N = 4) 
supersymmetry [Q], and for a particular gluon helicity configuration in pure Yang-Mills theory |15| ]. 
Subsequently, complete calculations have been performed for Bhabha scattering 16], general 2^2 
parton scattering in QCD [|l7], 18, the di-photon background to Higgs production at the LHC |2Q| , 



and light-by- light scattering [21]. 

Important technical breakthroughs, which allowed the more general calculations to proceed, 
included the reduction of the two loop momentum integrals appearing in the all-massless 2^2 
processes to a set of master integrals, and the evaluation of those master integrals as a Laurent 
series in e [22, 23, 24, 25, 26, 27]. Even more recently, the corresponding integrals where one of 
the four external legs is massive have been evaluated |28|], enabling the computation of the two-loop 
amplitudes for e+e~ annihilation into three partons [|2S|]. 

In our explicit verification of supersymmetry identities, we investigate the gg — > gg helicity am- 
plitudes 

^4(1^,2+3+4+), (1.1) 

where the subscript labels the particle species (g for gluon), and the superscript denotes the sign 
of the helicity. We use an "all-outgoing" helicity convention: if a given leg is incoming, then the 
actual helicity is the opposite of the superscript label. We study the helicity amplitudes ( |1.1| ) be- 
cause they are relatively simple (for two-loop amplitudes), and because supersymmetry is especially 
constraining. For both the ++++ and — +++ helicity configurations, we have computed all the 
gluon and fermion loop contributions that appear in QCD. For the ++++ configuration, we have 
also included scalars with both gauge and Yukawa interactions, to allow for a more extensive test of 
supersymmetry identities. 

It is important to be able to convert helicity amplitudes computed in the FDH scheme to the 
more standard CDR and HV schemes, and vice versa. A given dimensional regularization scheme 
has implications for regularization of both ultraviolet and infrared singularities. Both the ultraviolet 



and infrared aspects of scheme conversion have been extensively discussed at one loop [||, 11, 30], and 



to some degree at two loops [19]. The infrared aspects are not yet understood for arbitrary processes, 
as only the gg — > gg process was studied in ref. |i9| . Conversion from one scheme to another in the 
ultraviolet, though, is a process-independent procedure. It just amounts to relating the two different 
renormalized couplings implied by the two schemes, to a sufficiently high accuracy in perturbation 
theory. The FDH and DR schemes behave the same in the ultraviolet; their couplings are identical. 
The relation between the DR (or FDH) and MS versions of a s in QCD has been known to one-loop 
accuracy for some time |8], [31], [32], [llj . Here we extend the relation to two- loop accuracy. 

The first two coefficients of the QCD f3 function are scheme-independent (for analytic redefinitions 
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of the coupling), but the three-loop coefficient, &2> depends on the scheme. However, knowing the 
value of 62 hi the MS scheme and the two-loop relation between couplings, we can easily obtain 
the value of 62 in the DR (or FDH) scheme. 

This paper is organized as follows. In sect. [2] we review the supersymmetry Ward identities. We 
then present the rules for the FDH scheme in sect. [||. In performing the explicit two- loop calcula- 
tions, we did not use Feynman diagrams directly. Instead we computed (generalized) unitarity cuts 
of the amplitudes in all channels and to all order in e. Then the amplitude was reconstructed from 
the cuts. This cutting method is described in sect. ||. In sect. || we perform a color decomposition of 
the two-loop amplitude, into color structures multiplied by "primitive" functions which have been 
stripped of color. This decomposition makes it convenient to obtain either QCD or supersymmetric 
amplitudes. The primitive amplitudes used in our explicit check of the two-loop supersymmetry 
Ward identities are presented in sect. |6|. We verify the identities in sect. [7|. In sect. |8] we discuss ul- 
traviolet renormalization, the two-loop shifts in the couplings and amplitudes, and their implications 
for the three-loop (3 function. 



2 Supersymmetry Ward identities 



Helicity amplitudes in supersymmetric theories are subject to a set of stringent conditions imposed 
by the super-algebra: the S'-matrix supersymmetry Ward identities (SWI) |FJ, 13]. These identities 
allow us to distill the information contained in the supersymmetry algebra and apply it directly to on- 
shell S'-matrix elements. They hold in any supersymmetric theory. They also lead to relations among 
different components of amplitudes in any non-super symmetric theory, such as QCD, for which a 
re-adjustment of color representations and/or multiplicities makes the theory supersymmetric. At 
tree level and at one loop, supersymmetry Ward identities have been applied to QCD, either as 
checks or as computational aids [13, 34, 11, 35, ||. As we discuss in this paper, the same ideas can 
be applied at two loops. 



2.1 Derivation 

The derivation of the supersymmetry Ward identities from the super-algebra has been discussed in 
the literature in a number of articles and reviews [^, Q, so we describe it only briefly. Since we are 
interested in applications to non-supersymmetric theories we phrase the supersymmetry identities 
in terms of the component fields. The N = 1 super-algebra describing the action of the super-charge 
on the component gluon field g and gluino field A comprising the vector supermultiplet is 

[Q(p),g ± (k)]=Tr ± (k,p)\ ± (k), [Q(p),X ± (k)]=T^(k,p)g ± (k), (2.1) 

where k is the light-like momentum carried by the field, Q(p) is the super-charge contracted with a 
spinor for the (arbitrary) light-like vector p, and 

r + (k,p) = e[ P k] , r-(k, P ) = 9( P k) . 

Because T~ is proportional to a Grassmann variable 9 and T + is proportional to 8, the coefficients 
T + and T~ in a supersymmetry Ward identity are independent. We use the notation (k~\k^) = (ij) 
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and (kf\k- ) = [ij], where \kf) are massless Weyl spinors with momentum fej, labeled with the sign 
of the helicity and normalized by (ij) [j i] = = 2ki ■ kj. Note that (ii) = [ii] = 0. 

Since the super-charge Q(p) annihilates the vacuum, one can construct a typical N = 1 super- 
symmetry Ward identity in the following way: 

0= {0\[Q,gfgt\tgi}\0) =Tr ± (k 1 ,p)Af JSY (lf^,3i, 4+)-r + fe,pM 4 SUSY (1^2+ 3+,4+) 

-T'(k 3 ,p)Af JSY (l^,2+,3+, 4+) + r+(fc 4 ,p)A SUSY (1^2+ 3+4+), 

(2.2) 

where ^4.4 USY is a four-point amplitude in a supersymmetric theory, the integers refer to the leg 
labels, and the subscripts g and A to the particle species of the specified leg. Since the coefficients 
of r + and r~ are independent, the sum of the terms with T~ prefactors must vanish independently. 
(The r + terms in eq. (2.2) also vanish using gluino helicity conservation.) By choosing p = k% we 
obtain the identity, 

>tf USY (l± 2+ 3+ 4+)=0, (2.3) 

which is the main SWI that we will investigate at two loops in this paper. 

Using the super-algebra, one may systematically derive other identities. The nonvanishing am- 
plitudes with external gluons only are related to amplitudes containing external gluinos, e.g., 

4 SUSY,-,- 9 - o+ d +\ _ ( 13 ) .SUSY/,- 9 - o+ a + \ 
•^4 \ L g ! L \ : °A ' ^g I ~ /i o\ ^ \ y g i L g i °g > J > 



<12> (2.4) 



.SUSY n - o+ o- _ ( 24 ) /.SUSY/-,- 2 + o- 4 4 

4 V-'-A ' A ' A ' A / ~ (13) 4 ^ 9 ' Z 9 ' 6 9 ' 4 9 



The amplitudes on the right-hand side of eq. ( |2.4| ) have recently been computed at two loops . 
The amplitudes on the left-hand side are related to the QCD amplitudes for qq — > gg, qg — > qg, 
qq — ► qq, and qq — > qq, but these amplitudes have not yet been computed at two loops in the helicity 
formalism,Q and we shall not do so here. 



2.2 Lagrangian 

To allow us to separate out different supersymmetric combinations of amplitudes, we consider N = 1 
supersymmetric SU(N C ) gauge theory with two different types of matter content: 

1. rif identical chiral matter multiplets Qi, i = l,...,rif, transforming in the fundamental N c 
representation of the gauge group, and their N c partners Qi, with vanishing superpotential, 
W = 0. 

2. A single chiral matter multiplet $ = cj> a T a transforming in the adjoint representation, with 
superpotential W = ^g^Tr^ 3 , where Tr is an SU(N C ) trace. For convenience, we write the 
Yukawa coupling as and take £ to be independent of the gauge coupling g. 

1 They have been computed, and the SWI have been verified, at one loop [ pJ| . At two loops, the interference with 
the tree amplitude, summed over all external colors and helicities, has been computed JL7| |, but the conversion to a 
supersymmetric amplitude has not yet been performed. 
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Our SU(N C ) generators are normalized by Tr(T a T b ) = 5 ab . Because of the gauge symmetry, gluons 
couple to all other particles via the standard gauge theory interaction for fermions or scalars, with 
gauge coupling g. The couplings of gluinos or matter fermions may be conveniently extracted from 
the component expansion of the interaction Lagrangian . 

The interaction terms in the case of fundamental matter are 

n f 

l4 und = igJ2U A * Ta ^ - i>iT a A*)\ a - A°(4T°& - TpiT a Ai) 

(2-5) 

g 2 , «r 



^h(A*T a Ai — AiT a A*)\ , 
i=i ' 



4 

following the notation of Wess and Bagger [|37j , chapter 7. Here A = X a T a is the gluino; Ai and ipi 
are, respectively, the scalar and fermionic components of Qi, while A, and ipi are the corresponding 
components of Qi. 

In the case of adjoint matter, the interaction terms proportional to the gauge coupling are, 



L adjoint = _ ig TT{[A * jm+l g Tv{[A,#]\) - ^ £ M[A A*}T°) 



4 

-ig Tr([A*,ip]X) + ig Tr([A, %b}\) - — Tr (AA*AA* - A 1 A* 



(2.6) 



The terms arising from the superpotential W are 

L adjoint = _^ Tv ^ A _ g £* Tv ^ A * _ 5 2|^|2 ^ Tr ^ T ^ Tt(A* 2 T° 



-gi TrWA-g? Tr^A* ~ 9 2 \C\ 2 



TiA 2 A* 2 - — TtA 2 TtA* 2 



N. 



c 



(2.7) 



2.3 Amplitude decomposition by particle content 

Supersymmetry can be exploited in non-supersymmetric theories such as QCD [jl3j, |||, ^] by observing 
that appropriate linear combinations of quantities occurring in QCD amplitudes are in fact super- 
symmetric. As an especially simple example, at tree level an n-gluon amplitude is automatically 



supersymmetric [13]. Since the fermions do not appear in intermediate states, the fermions in the 
theory might as well be in the adjoint representation, i.e., the theory might as well be pure N = 1 
super- Yang-Mills theory. Thus, by virtue of eq. (|2.3D, the four-gluon tree amplitude in QCD satisfies, 



^QCD tree (1 ± 2 + 3 + 4+) = . (2.8) 

At loop level the application of supersymmetry in QCD is clearly more intricate, because all 
particles in the theory can circulate in the loops. Nevertheless, one can still use supersymmetry to 
relate different contributions. 

At one loop, the four-gluon amplitude in supersymmetric QCD with rif matter multiplets — case 



1 in sect. 2.2 — has the structure, 

Af =1 = AT ctor + n f Af atteI , (2.9) 
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where the contribution ^4| ector is that of a vector multiplet consisting of a gluon and gluino, and 
^matter j g fag^ Q f a matter multiplet consisting of a quark and squark. Identities may be obtained 

corresponding to 



from each term in this decomposition. For example, by setting nj = in eq. (|2. 
pure N = 1 super- Yang-Mills theory, from the SWI (2.3) we have, 



A 



gluon loop / -I ± 



I 31 2 



3 + 



--4: 



gluino loop/-! ± 



I 31 2 



3 + 



(2.10) 



which relates the gluon loop contribution depicted in fig. ||(a) to the fermion loop contribution 
depicted in fig. |l](b). To apply this identity to QCD, one uses group theory to relate the contribution 
of an adjoint representation fermion (gluino) in the loop to that of a fundamental representation 



fermion (quark), as we shall describe further in section 5.1. 

Similarly, by considering the contribution from an N = 1 matter multiplet, consisting of a fermion 
and a scalar transforming in the same representation of the gauge group, i.e., the nj-dependent term 
in eq. (|2.9Q , we obtain 



^ caJarloop (l± 2+ 3+ 4+) 



\ fermion loop / -i ± n+ 
4 \ L g i Z g 



(2.11) 



which relates the fermion loop contribution in fig. [j](b) to the scalar loop contribution in fig. [j](c). 
By supersymmetry, the number of physical states in the scalar loop matches that in the fermion 
loop. For a scalar in the adjoint representation of SU (N c ), paired with a Majorana gluino (in N = 2 
supersymmetry, say), this amounts to 2 x — 1) states. For a fundamental representation scalar, 
paired with a quark (a Dirac fermion), it amounts to 4 x N c states. (At the one- loop level, considering 
adjoint matter with a superpotential — case 2 in sect. |2.2| — leads to no new identities from the 
four-gluon amplitude.) 




(b) 



(c) 



Figure 1: Sample diagrams contributing to: (a) a gluon circulating in the loop, (b) a fermion circulating in 
the loop, and (c) a scalar circulating in the loop. 



These considerations extend to two loops. In the case of rif matter multiplets in the fundamental 
representation — case 1 in sect. |2.2| — the four-gluon amplitude takes the form, 



A N=1, fund = ^vector + nfA ™*<*W + ^matter^J _ (2 12) 



In the case of one adjoint matter multiplet with a superpotential — case 2 in sect. 2.2 — the 
four-gluon amplitude takes the form, 

^A r =ladj _ ^vector _|_ ^matter(l) ^matter(2) |£|2^jYukawa ^ ]^ 
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The quantities ^4 matter M are "dressed" differently with color in the two cases, as will be explained 
in sect. |5[ We have subdivided the contributions according to their dependence on the number 
of matter multiplets and on the couplings. Identities may be obtained from each term in the 
decomposition ( 2.12 ), and from the ^-dependent term in eq. fl2.13|) , because nf and £ are independent 
parameters. Representative diagrams contributing to each of the four independent supersymmetric 



components A\ ector ) A 



vector ^ matter (1) .matter (2) 



A 



and ^Yukawa are depicted in figs. @(a)-(d). 




4 1 



4 1 





(a) 

3 2^ :! 2 



1 



1 



(b) 







(c) 



2 \ 




(d) 

Figure 2: Representative diagrams contributing to the supersymmetric amplitudes (a) A 



vector ^matter (1) 

(c) ^™ attcr ( 2 ) 5 anc j (d) A% ukawa in eqs. ( 2.12Q and (2.13). The curly, solid without arrows, solid with arrows, 



and dotted lines represent gluons, gluinos, quarks, and scalars, respectively. 



For the helicity configurations ++++ and — +++, four independent identities are obtained by 
setting eqs. ( 2.12| ) and (2.13) to vanish using the SWI fl2.3p . Each such identity actually generates 
several equations, once the amplitudes are further decomposed according to the colors of the external 
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gluons. The equations stemming from .A vector (=t, +, +, +) = 0, which only involve gluons and gluinos, 
can be applied to QCD once one understands the group theory relations at two loops between the 
contribution of fundamental and adjoint representation fermions. The remaining identities have 
scalar particles propagating in the loops as well. A systematic discussion of all these relations, based 
on "primitive" or "color-stripped" amplitudes, will be presented in sections |5|]7| 

3 Supersymmetric regular izat ion 

For supersymmetry identities to hold directly, a necessary condition is that the regularization should 
not alter the number of bosonic states relative to the number of fermionic states. The conventional 
and 't Hooft-Veltman variants of dimensional regularization [Q, are incompatible with supersym- 
metry precisely because they alter the balance of bosonic and fermionic degrees of freedom. In the 
four-dimensional limit, manifest supersymmetry will not generally be recovered in these schemes due 
to divergences in the amplitudes. An 0(e) discrepancy between the number of bosonic and fermion 
states can be multiplied by an 1/e singularity, leaving a supersymmetry- violating remainder even 
as e — > 0. If the only divergences are ultraviolet in nature, it is possible to repair such violations 
by adding suitable finite counterterms, order by order in perturbation theory (see e.g. ref. |32|j). 
However, it is clearly desirable to avoid this situation if possible, particularly in theories with a large 
number of coupling constants. Also, we are interested in on-shell scattering amplitudes for theories 
with severe infrared divergences, and here local counterterms will not suffice. 

3.1 Dimensional reduction scheme 

The most widely used scheme for preserving supersymmetry is the dimensional reduction (DR) 
variant of dimensional regularization.^] The rules for dimensional reduction follow from viewing it 
as a compactification of a four-dimensional theory to D < 4. The rules for dimensional reduction 
are ||: 

• As in ordinary dimensional regularization, all momentum integrals are integrated over D- 
component momenta. Any Kronecker resulting from the integration are D-dimensional. 
(This is necessary for maintaining gauge invariance.) 

• All indices on the fields, and on corresponding matrices coming from the action, are treated 
as four-dimensional indices. 

• Since D < 4 always, any four-dimensional Kronecker S^ u contracted with a D-dimension 
momentum pi D ^ yields a D-dimensional momentum, Sjx pu = pju • Similarly, for any four- 
dimensional vector, , the dot product with a D-dimensional vector yields a D-dimensional 

dot product, £^pjP' 1 = ej^pjL . In general, dot products of four-dimensional vectors with 
2e-dimensional ones will not vanish. 

2 As is customary in essentially all variants of dimensional regularization, we treat the fermions as four-dimensional 
by letting the Dirac trace of the identity be tr(l) = 4. 
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The first rule is necessary for preserving gauge invariance and the second for preserving super- 
symmetry. It is the third rule which defines the regularization as dimensional reduction. In the DR 
scheme a four-dimensional vector may be viewed as a combination of the D = 4 — 2e vector plus 
a set of 2e scalars. In non-supersymmetric theories in the DR scheme it is especially important to 
keep track of the distinction between the vectors and the 2e scalars because of their differing renor- 
malization properties. Moreover, in the non-supersymmetric case it is essential to keep evanescent 
couplings and operators ||. 

The rule &&Pp = £ ^d)P^ 1S a ^ so awkward to handle in the presence of explicit four-dimensional 
polarization vectors e^. Such vectors are encountered when evaluating helicity amplitudes in (for 
example) the spinor helicity formalism Q. Intuitively, for D < 4, there are less than 2 spatial 
directions transverse to the gluon direction, so one cannot really perform a rotation in this transverse 
plane, as required to define a helicity eigenstate. 

3.2 The four-dimensional helicity scheme 

Motivated by the desire of having a supersymmetric scheme whose rules are more closely related to 
the more conventional 't Hooft-Veltman scheme, and more compatible with the helicity formalism, 
we define the four- dimensional helicity scheme, which has already been used in a number of one-loop 
calculations || HH> H> ■ ^he essential difference between the DR scheme and the FDH scheme is 
that in the former case the rules for dot products follow from taking D < 4 while in the latter one 
they follow from taking D > 4. 

Our rules for extending the four-dimensional helicity scheme to two loops are as follows: 

• As in ordinary dimensional regularization, all momentum integrals are integrated over D- 
component momenta. Any Kronecker 5"'s resulting from the integration are D-dimensional. 
(This is necessary for maintaining gauge invariance.) 

• All "observed" external states are left in four dimensions; their momenta are also four- 
dimensional. (In QCD the "observed" states refer to the external states appearing in the hard 
part of the process described by Feynman diagrams, ignoring any subsequent hadronization.) 
Because D > 4, we may view this rule as choosing momenta and polarizations to lie solely 
in a four-dimensional subspace. In this way it is natural to use helicity states for "observed" 
particles. 

• All "unobserved" internal states are treated as D s dimensional, where D s > D in all intermedi- 
ate steps. The "unobserved" states include virtual states in loops, virtual intermediate states 
in trees (which may be attached to loops), as well as any external states which are in collinear 
or soft parts of phase space. Any explicit factors of dimension arising from the Lorentz and 
7-matrix algebra should be labeled as D s , and should be kept distinct from the dimension D 
describing the number of components of the loop momenta. 

• Since D > 4, for any four-dimensional vector (such as an "observed" polarization vector or 
momentum), e^, the dot product with a L>-dimensional vector yields a four-dimensional dot 
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product, £(4)?ty = £(4)Pft • I* 1 general, dot products of four-dimensional vectors with (— 2e)- 
dimensional ones always vanish. 

Though the rules for various dot products are constructed with D > 4 in mind, at the end the 
expressions are analytic functions of D and D s which can be continued to any desired region. By 
setting D s = D these rules are precisely the ones for the 't Hooft-Veltman scheme. The FDH 
scheme is specified by taking the parameter D s — > 4, after all Lorentz and 7-matrix algebra has been 
performed. In performing the 7-matrix algebra, the 't Hooft-Veltman prescription [2] for 75 should 
be used (75 commutes with 7^ when the index \x lies outside of four dimensions). 

A feature of the FDH scheme that makes it useful in QCD amplitude computations is its simple 
relation to the widely used 't Hooft-Veltman scheme. By keeping track of the D s parameter when 
performing calculations one can easily switch between the 't Hooft-Veltman and the FDH scheme. 
In section ||, we will quote results in both the FDH and the 't Hooft-Veltman scheme, leaving D s as 
a free parameter. In fact, one can define a continuous class of schemes by setting D s = 4 — 2c5r, as 
in eq. ( |8.5[ ) below, though we see no particular utility to schemes other than FDH (5r = 0) and HV 
(*fl = l). 

For this approach to be sensible, the coefficients of each power of D s must be separately gauge 
invariant. In any theory with either gauge or Yukawa interactions, the terms with D s in a given 
Feynman diagram can be mapped to another diagram where (fictitious) scalar lines replace some of 
the gluon lines in the diagram. After summing over diagrams, the terms with D s to a certain power 
are proportional to a sum of diagrams containing a certain number of fictitious scalar loops; the sum 
is gauge invariant because it corresponds to a physical amplitude. 

At one loop, the supersymmetry preservation properties of the FDH scheme have been verified 



in a number of papers [11, One of the aims of this paper is to provide explicit examples 

demonstrating that the FDH variant of dimensional regularization preserves the SWI fl2.3|) through 
at least two loops. 

4 Cutting method 

As mentioned in the introduction, we did not use Feynman diagrams directly to compute the explicit 



two-loop helicity amplitudes. Instead we used a cutting method |35], |l, ^j, which has been applied 
previously to a number of one-loop calculations, including the corrections to e + e~ — > 4 partons |4C| [, 
and more theoretical studies, such as the construction of infinite sequences of maximally helicity 
violating amplitudes fi^] and the investigation of the divergence structure of supergravity j4jj. 
More recently, it has also been used to produce two-loop 2 — > 2 scattering amplitudes in super- Yang- 
Mills theory, QCD and QED [14, 15, 21, [l9|]. In sect. || we will present the four-gluon amplitude 
for the ++++ helicity configuration, in both QCD and supersymmetric theories, obtained via the 
cutting method. We use this amplitude to investigate the supersymmetry identities and associated 
regularization issues. We have also computed, by the same techniques, the — h++ amplitude com- 
ponents which do not involve scalars; they too satisfy the supersymmetry identities. The cutting 
method can help clarify the unitarity and gauge invariance of the regularization procedure, because 
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the basic building blocks for loop amplitudes are gauge invariant tree-level S'-matrix elements. 

The cutting method amounts to an extension of traditional unitarity methods [44]. Traditional 



applications of unitarity in four dimensions, via dispersion relations, often suffer from subtraction 
ambiguities. These ambiguities are related to the appearance of rational functions with vanishing 
imaginary parts, R(Si), where Si = {s,t, u, . . .} are the kinematic variables for the amplitude. 
However, dimensionally-regulated amplitudes for massless particles, as we consider here, necessarily 
acquire a factor of (— Si)~ e for each loop, from the loop integration measure / d 4 ~ 2e L and dimensional 
analysis. For small e, we expand (— 5j) _<E R(Si) = R(Si) — e ln(— Si) R(Si) + 0(e 2 ), so every term 
has an imaginary part (for some Si > 0), though not necessarily in those terms which survive as 
e — > 0. Thus, the unitarity cuts evaluated to 0(e) provide sufficient information for the complete 
reconstruction of an amplitude through 0(e ), subject only to the usual prescription dependence 
associated with renormalization. The subtraction ambiguities that arise in traditional dispersion 
relations are related to the non-convergence of dispersion integrals. A dimensional regulator makes 
such integrals well-defined and correspondingly eliminates the subtraction ambiguities. In a sense, 
we use dimensional regularization as a calculational tool, beyond its usual role as an infrared and 
ultraviolet regulator. 

It is useful to view the unitarity-based technique as an alternate way of evaluating sets of ordinary 
Feynman diagrams. It does this by collecting together gauge-invariant sets of terms which correspond 
to the residues of poles in the integrands. The poles are those of the propagators of the cut lines. 
This corresponds to a region of loop-momentum integration where the cut loop momenta go on shell 
and the corresponding internal lines become the intermediate states in a unitarity relation. From 
this point of view, we may consider even more restricted regions of loop momentum integration, 
where additional internal lines go on shell (and, if they are gluons, become transverse as well). This 
amounts to imposing cut conditions on additional internal lines. 

Besides the more traditional two- and three-particle cuts one can define "double" two-particle 



generalized cuts [15] for a two-loop four-point amplitude. An example of this quantity is illustrated 



in fig. 0(a), and written in terms of on-shell tree amplitudes as, 



A^ 



2-loop 



, , ♦= E ^ ee (^,-e 2 ,-£i)xA t i ee (e l ,e 2 ,-e 3 ,-e 4 )xA t i ee (e 4 ,£ 3 ,3,A), (4.1) 

2x2-cut ' — ' 

physical 



where the on-shell conditions if = are imposed on the £4, i = 1,2,3,4 appearing on the right- 
hand side. This equation should not be interpreted as trying to take "the imaginary part of an 
imaginary part" . Rather it should be understood in the sense of the previous paragraph as supplying 
information about the integrand of the two-loop amplitude. It supplies only part of the information 
contained in the usual two-particle cut, which effectively imposes only two kinematic constraints 
on the intermediate lines. However, it is simpler to evaluate because it is composed only of tree 
amplitudes. There are, of course, other ways to cut the two-loop amplitude to obtain trees. For 
example, in fig. 0(b) a different arrangement of the cut trees is shown, 



A 4 



2-loop 



hv-cut 

physical 
states 



J2 AT»0-,ei,e 2 , -h) xAT e (2, -e^-e^-ej x < ee (^ 3 ,3,4). (4.2) 
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The combined set of double two-particle cuts provides all information present in the ordinary two- 
particle cuts (where a single pair of lines is cut, and a loop amplitude is still present), thus obviating 
the need to evaluate such cuts. 





Figure 3: Two examples of s-channel double two-particle cuts of a two-loop amplitude, which separate it 
into a product of three tree amplitudes. The dashed lines represent the generalized cuts. 



The full amplitude, including all color factors, may be obtained by combining a suitable set of 
generalized cuts into a single expression whose cuts match the explicitly calculated cuts. At two 
loops, it is sufficient to evaluate all the double two-particle cuts with the topologies shown in fig. ||, 
plus the "standard" three-particle cut, shown in fig. and given by 



A 



2-loop 



3-cut 



(4.3) 



physical 
states 



For the identical-helicity amplitudes discussed in this paper, the integrands are sufficiently simple 
that we can combine all double two-particle cuts and three-particle cuts into compact integrands 
containing no cut restrictions. We present these compact integrands in sect. 0. 




Figure 4: The standard three-particle cut of a two-loop amplitude. 



To obtain and verify the compact representations of the amplitudes presented in this paper, 
it proved useful to compare numerically, at a number of random kinematic points, two different 
representations of the cut integrands ("raw" and simplified), before performing any loop integrations. 
This comparison is only simple to implement when the number of dimensions D is an integer. Integer 
values for D may seem at odds with dimensional regularization, which requires expressions to be 
evaluated in non-integer numbers of dimensions in order to analytically continue to D = 4. However, 
a cut integrand contains no explicit dependence on the dimension D, only that implicit in the dot 
products of loop momentum vectors. It is therefore sufficient to verify the cut integrands numerically 
for integer values of D. In so doing, we keep D s an analytic parameter, independent of D; i.e., we 
verify the cut integrands for all values of D s . One should take care that the number of dimensions 
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is no smaller than the number of independent vectors in the problem. At two loops, there are two 
independent loop momenta with non-vanishing components in the extra dimensions. Hence the 
extra-dimensional subspace must be at least two-dimensional, for a minimum total dimension of 
D = 6. If two cut integrands agree numerically in six or more dimensions (with D s left arbitrary), 
then the lack of any explicit D dependence ensures that they are identical for any value of D. 

5 Color decomposition of QCD amplitudes 

Gauge theory scattering amplitudes have a rich color structure. A number of different color de- 
compositions have been used to organize this structure, particularly at the tree and one-loop 
level |45| , ^ |4(], H ^7]]. In general, an amplitude is color decomposed by writing it as a sum of 
terms. Each term is the product of a "color structure" and a "primitive" (or "color-stripped") 
amplitude. A color structure is a tensor in color space, but is independent of the momenta and 
polarizations of the external states. A primitive amplitude, on the other hand, contains no color 
indices or group theory information; it is a function only of the kinematic variables. 

For gluon scattering amplitudes in SU(N C ) gauge theory, "trace-based" decompositions have 
been used frequently fig , |4], g]. Here the color structures have the form TV(T ai i ■■■T ain ), 
TV(T a4 i ...T aim ) r Ti(T a%m + 1 ...T ain ), etc., where T a is a generator in the fundamental represen- 
tation. However, a color decomposition based on the Lie algebra structure constants f abc [^7]] is 
more suitable for our purposes. Note that the set of color structures may be complete, that is, it 
may form a linearly independent basis; or it may be overcomplete, with its elements obeying linear 
relations. (Usually the completeness is defined with respect to an arbitrarily large value of N c ; for 
small N c it may degenerate.) If the set is complete, then the primitive amplitudes are uniquely 
defined; if it is overcomplete, then there is some freedom in defining them, although there may still 
be a natural, symmetric way to do it. 

The utility of primitive amplitudes in the context of supersymmetry is that they can serve 
as building blocks for amplitudes in disparate theories, where the matter transforms in different 
representations of the gauge group. If one of the theories is supersymmetric, one can express super- 
symmetry relations in terms of primitive amplitudes. Then the supersymmetric properties of the 
primitive amplitudes can be applied to non-supersymmetric theories such as QCD. 

In this section, we first review tree and one-loop color decompositions, and remind the reader how 
they allow supersymmetry Ward identities to be applied to non-supersymmetric theories. Next we 
proceed to two loops. We organize the color and kinematics of the two-loop gg — > gg amplitudes in 
a representation convenient for discussing the supersymmetry identity ([D|). For the case where all 
particles in the loops are in the adjoint representation, the two-loop color decomposition was given 
previously p3[ pl|] . Here we review this result and extend it to the case of matter in the fundamental 
representation. The sets of color structures used in both decompositions are overcomplete, and 
therefore the primitive amplitudes they define are not unique. However, the overcompleteness gives 
us the freedom to find very symmetric and compact forms for the primitive amplitudes, which do obey 
simple supersymmetry relations, at least for the special case of the +- 1— I— I- helicity configuration. 
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5.1 Color organization of tree and one- loop amplitudes 



In order to explain the utility of primitive amplitudes in applying supersymmetry identities to QCD, 
first consider tree amplitudes. We employ an / a6c -based (not a trace-based) color decomposition |47| ], 
because it most closely matches the one we shall use at two loops. Using the Jacobi identity for the 
structure constants, the four-gluon tree amplitude can be decomposed as a sum of two terms, 

AT e (l 9 ,2 g ,3 9 A 9 ) = 9 2 E (i ?aCT(2) i ?aCT(3) )a 1 a 4 < ee (l 9 ^(2 9 ),a(3 9 ),4 9 ), (5.1) 

(7G5 2 

where the primitive amplitudes A^ ee do not contain any color information. Here S2 is the set of 
two permutations of the set {2,3}, while (i ?a )fe c = if bac is an SU(N C ) generator for the adjoint 
representation, 

jabc = ^jabc = - i T [ ([T a ,T*]T e ) , (5.2) 

where f abc are the usual SU(N C ) structure constants, and T a are generators for the fundamental 
representation, normalized so that Tr(T a T b ) = 6 . Because the two color structures in eq. (5.1) 
are linearly independent, the SWI fl2.8|) for „4^ CD tree applies separately to each primitive amplitude 
in eq. (|5.1|) . 



Next consider the two-gluino two-gluon tree amplitude in a supersymmetric theory. Because 
gluinos are in the adjoint representation, this amplitude satisfies a color decomposition equivalent 
to eq. Q, 

AT (1a,2 9 ,3 9 ,4 a ) = g 2 £ (FM 2 ) F*°W) ailH A¥»(l x , a{2g ), a(3 g ), 4 A ) . (5.3) 

Compare this decomposition to the one for the two-quark two-gluon amplitude in QCD, 

AT e (l q ,2 9 ,3 9 ,4g) = 9 2 E {T a ^T a ^)^AT e {l q ,cj{2 g ),a{3 g )A,) , (5-4) 

<res 2 

where the quarks, and T ai , are in the fundamental representation. The crucial point is that the 
A^ ee appearing in eq. (|5l^) are in fact identical to the Alp 6 appearing in eq. (|[x4|), as can easily be 
verified using Feynman diagrams. (In the double-line formalism for SU(N C ), the color line running 
between the two gluinos is merely "stripped off" to obtain the quark amplitudes.) Thus, there is no 
real distinction between gluinos and quarks as far as the primitive amplitudes are concerned. Any 
SWI that holds for the primitive gluino amplitudes will hold for the quark ones. 

Similar color decompositions hold at one loop. Consider the four-gluon amplitude. When the 
particle circulating in the loop is in the adjoint representation, e.g. a gluon or gluino, the decompo- 



sition is [47 1 



^adjoint i°°P( 1)2 ,3,4) = <? 4 V Tr(F a -( 1 )F a -(2) j p a -(3)F a -( 4 ))^" loop (a(l),c7(2),cj(3),cj(4)) 

-a£S 4 /Z 4 /n 

(5.5) 

Here we have suppressed the gluon g labels on the external legs; a runs over the set of permutations 
of {1,2,3,4}, after removing those equivalent under a cyclic Z4 permutation or the reflection 1Z: 
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{1,2,3,4} — * {4,3,2,1}. In comparison, the contribution of a fundamental representation particle 
in the loop, e.g. a quark, is decomposed as 



4 und - loop (l,2,3,4) =g A 



]T Tr(TMD . . . TM4) ) A l - Xoop {a{l), <r(2), a(3), <r(4)) 

■CTGS4/Z4 



(5.6) 

Note that in both eqs. ( |5.5| ) and fl5.6j), the color factors may be read off of the appropriate one-loop 
"parent" diagram shown in fig. ||. One simply assigns a group theory factor of f abc or (T a )j for each 
vertex, and a factor of S ab or 5j for each internal line, then performs the index contractions. We 



shall use this diagrammatic representation of the color factors at two loops in sections |5.2| and |5.3| . 

Once again, using Feynman diagrams, it is straightforward to demonstrate that the primitive 
amplitude A\ loop appearing in the fundamental representation case (5.6) is exactly the same object 



appearing in the adjoint representation case ( |5.5|) , provided that the type of particle in the loop 
(fermion or scalar) is the same in both cases. 

Because of this identification, QCD quark loop primitive amplitudes obey SWI. For example, let 
us apply eq. ( |2.3| ) to pure N = 1 supersymmetric Yang-Mills theory (SYM). Linear independence of 
the color factors appearing in eq. ( |5.5[ ) implies that the primitive amplitudes satisfy 

4 YM 1 - loop (l ± ,2 + ,3 + ,4+) = Af uon loop (l ± ,2+,3 + ,4+) +4 ermion lo °P(i±, 2+, 3+ 4+) = 0, (5.7) 

where ^4| luon loop represents the gluon loop contribution to A\ loop in eq. (|5.5j), shown in fig. [TJa, and 

^fermion loop , , q 1 • 1 ±. -i ; ■ <l-loop i • n 171 r, ^fermion loop 

j4 4 represents the giumo loop contribution to A A , shown m hg. |l|b. Because A A 

is the same for an adjoint or a fundamental representation fermion (although it gets inserted into a 



different color decomposition formula, (5.5) or ( |5.6| ), in the two cases), the identity 



^giuon i°op( 1 ± )2 +,3+,4+) = -A^ mion loop (l ± ,2+,3+,4+) (5.8) 

holds even in QCD, where the fermion would be a quark. Continuing along these lines, consider the 
contribution of a chiral multiplet in super-QCD, consisting of a fermion and a scalar, to eq. ( |2.3| ). 
One obtains 

^scalar loop (1 ± ^ 2 + ^ 3+ ^ 4+) = _^fermion loop^ 3+ 3+ 4+) ^ ^ 



even for a non-super symmetric gauge theory. As discussed in section 2.3, one must take the number 
of scalar states to match the number of fermion states on the two sides of eq. (|5,9| ). 

It is not difficult to verify that the SWI $)J^) and ( |5.9[) are indeed satisfied. For the identical- 



helicity case (+++- 1-), the one-loop amplitudes have a compact representation [41, 48], similar to 
the one we shall use at two loops, 

Af luon loop (l+, 2+ 3+ 4+) = (D s - 2) ?rll- loop [\*} , 

^fermion loop (1 + ^ 2+ ^ 3+ ^ 4+) = _ 2 P_ ^l-loop ^ ^ (5 _ 1Q) 
^scalar loo P(1 + 2 +, 3 +, 4+) = 2 ^ l\' X °° V [A*] , 

corresponding to the three representative diagrams in fig. |. Here 



p - i' 12 '' 34 ' (5.11! 

H (12) (34) v ' 
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is a ubiquitous spinor product prefactor, which is a pure phase, and is totally symmetric under 
permutations of the external legs {1,2,3,4}. The one-loop integral appearing in the amplitudes is 

T i-ioo P r A 4i _ f d D p 

4 [ pl J (2tt)d p 2( p - ki) 2 (p -h- k 2 y{p + k A y (5 12) 

where we have split the loop momentum p = pm + X p into its four-dimensional components pw 
and its (— 2e)-dimensional components X p , with A 4 = (X p • Xp) 2 ■ (The value of x\ loop [Ap] is known 



through 0(e) [15].) Full amplitudes, including all color factors, are obtained by substituting the 



primitive amplitudes in eq. ( 5.10|) into eqs. ( |5.5|) and (|5 



From eq. (5.10) we see that the SWI ( |5,8D and ( |5.9| ) hold to all orders in e in the FDH scheme, 
where D s = 4. Note that in the HV scheme, where D s = 4 — 2e, the SWI do hold at one loop in the 
limit e — > for the ++- 1-+ (and — h++) helicity configurations 0; however, the identities for the 



1-+ and — H h amplitudes are known to be explicitly violated, even in this limit [11]. 

5.2 Purely adjoint representation two-loop color decomposition 

The two-loop color decomposition we use is similar to the tree and one-loop decompositions described 
above. For the case where all particles are in the adjoint representation, the complete amplitude is 
given by 

4 adjoint _ 6 [^D; ^D; , ^Dj *T>i , n(cyy A \ 

•^X —9 2^ [°1234 ^X1234 L/ 3421 Ji X3A21 H^vj • (5 13) 

{DO 

Here C Di are color factors, while A Di are primitive amplitudes. The notation "+ C(234)" instructs 
one to add the two non-trivial cyclic permutations of {2, 3, 4}. To simplify the formulae, we adopt a 
diagrammatic representation of the color factors, in terms of a set of "parent" diagrams, depicted in 
figs. |^-||. The label X = G, S,F,M refers to whether the contribution is respectively from gluons, 
adjoint scalars, adjoint fermions, or mixed fermion-scalar loop contributions, corresponding to figs. ||, 
|B], ^ and ||, respectively. Which of these contributions appear in the full amplitude depends, of course, 
on the matter content of the theory under consideration, through equations like ( 2.12j ) and ( |2.13| ). 



Each color factor C D ' in eq. (5.13) is specified by a parent diagram, where the corresponding 



label Dj is shown in parentheses in the figure. For example, the fermion loop contribution color 
factor 01234 is found from fig. §(Pi). The color factor is calculated from the parent diagram by 
associating each vertex and internal line with a color tensor, and then performing the internal index 
contractions. For the case where all particles are in the adjoint representation, the internal lines are 
"dressed" with factors of 5 ab . The vertex dressing rules are as follows: 

• A gauge 3-point vertex is dressed with if abc = Tr([T a , T b ]T c ). 

• A superpotential 3-point vertex is dressed with d abc = Tr({T a , T b }T c ). 

• The 4-point vertex in the double-scalar loop diagram P5 in fig. ^ is dressed with 

* ±i 2 (f abe f ecd + pcejeda-j j n ^ e case Q f a gauge coupling (D term), and 
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* d bce d eda in the case of a superpotential coupling (F term). 



The structure of the 4-point vertices follows from the Lagrangian terms ([^(j) and ( j2.T|) , respectively. 
Color factors with the legs ordered differently from the figures, such as Cj^, are obtained by 
appropriate relabeling of the external legs. In the fundamental representation case to be discussed 
in section ^1^, the same color factor rules will hold, after correcting for the different representations; 
for example, b\ should clearly be used for a fundamental line, and (T a )j for a three-point vertex with 
two fundamental lines emanating from it. 

The primitive amplitudes A Di are defined as the coefficients of the color factors in eq. ( |5.13| ). 
Actually, in contrast to the tree and one-loop cases, this prescription does not completely specify the 
primitive amplitudes, because the color factors in eq. ( |5.13| ) are not linearly independent (see below). 
Also, the quantities A^ ee and A\ loop can be given gauge-invariant definitions, as sums of complete 
sets of color-ordered Feynman diagrams. Such a definition fails here. For example, the pure glue 
contributions ^4^1234 anc ^ ^G*i432 both contribute at leading order in N c to the cyclic color ordering 
1234, and thus only the sum of them has a gauge-invariant definition. In any case, the primitive 
amplitudes A Di that we present in sect. || are arranged to fulfill eq. ( 5.15j ). In addition, each such 
A Di does include the corresponding parent diagram from figs. in its definition; it also includes 
pieces of other "daughter" Feynman diagrams (not shown in the figures) , which are typically shared 
with other primitive amplitudes. 





1 (P) 4 1 (NP) 4 

Figure 5: Parent graphs for the pure gluon G contributions. 



For the purely adjoint color representation, in the gauge coupling case for which all vertices are 
built out of factors of f abc , an equivalent but much simpler color decomposition has been given 43, 
|i"5| ]. This decomposition is in terms of just the color factors P and NP shown in fig. |5[ Subdividing 
the parent diagrams Dj into planar ones Pj and non-planar ones NPj, we have 

.adjoint e/'spLoP a p 1 , r p A Pi 1 , v Tr NP a np * +r NP 4 NP * 1 4- rd^X 

—9 \ 2^ |_°1234 yl X1234 ' ^3421 /i X342l\ + [ U 1234 /1 X1234 + U 3421 /1 X342lJ + U ^«J4| > . 

L {P ! } {NPJ ' } 

(5.14) 

In the pure gluon indicated in fig. |5|, we lump all of the Pj primitive amplitudes together, 

and similarly for the NPj ones, so that the actual decomposition is 

•A-Q* = 5 f6 |cf234 ^G1234 + ^3421 ^4g*3421 + ^1234 ^1234 + ^3421 ^G3421 + C(234)| . (5.15) 

The equivalence of eqs. ( 5.13| ) and ( |5.14| ) is not completely manifest. While it is clear that 

Ly 1234 ~~ °1234 ) u 1234 — ° 1234 ) 1 — 1,4,0, {O.IU) 
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it is also true that 

^1234 1- ^1234 T ^1234 • 

Nevertheless, the symmetry properties of the primitive amplitudes presented in sect. [|| are such that, 



u 1234 T u 1234 T u l 



after summing over the permutations in eqs. ( 5.13| ) and (|5.14 ), the two forms are equivalent. The 



P4 primitive amplitudes, after loop integration, all turn out to be proportional to the same function, 
which is antisymmetric under exchange of legs 1 and 2 (or 3 and 4) . In the permutation sum ( |5.14| ) , 
this property induces an antisymmetric projection on the double planar box (P) color factor, which 
removes certain unwanted (subleading-color) terms, and renders it equivalent to the double triangle 
(P4) color factor ||15|| . Similarly, the double-scalar P5 primitive amplitude is symmetric under ex- 
change of legs 1 and 2 (or 3 and 4), and the symmetric projection on its C 12 5 3 4 color factor in eq. ( |5.13 ) 



renders it equivalent to C p 2 34- We emphasize that the simplified color decomposition (5.14) applies 



only to the pure gauge coupling, adjoint representation case, and not when either Yukawa couplings 
or fundamental representations are present. 

The symmetries of the color factors C p and C NP can be read off the diagrams. The ones required 
in the sum (|5.14|) are 

r<P — r*P r<NP _ ^np /c -1 7 \ 

u 4321 — Ly 1234 1 1 243 — u 1234 • \°- L <) 

The corresponding planar and non-planar primitive amplitudes share the same symmetries with their 
associated color factors. Due to these symmetries, .A^ omt (l + , 2 + , 3 + , 4 + ) has the required total (64) 
permutation symmetry, even though only six permutations appear explicitly in eq. ( jug) . Although 
diagrams with differing color factors contribute to each primitive amplitude, roughly speaking, gauge 
invariance dictates that in the final gauge invariant expression, all diagrams follow the lead of the 
parent diagrams. The final permutation sum then ensures that diagrammatic contributions with 
seemingly incorrect color factors receive the correct ones after assembly. 

The color representation in eq. ( |5.14j ) is not unique, although it is a particularly symmetric one. 
Indeed, the twelve color factors appearing in eq. ( 5.14] ) or ( p. 15 ) satisfy a set of seven linear relations, 



p 

1342 


r P 
°3421 


^NP ,-rNP 
— °1342 ~~ °1234 i 


n P n P 
°1423 — °4231 


r^p 

— u 1423 


^NP 
~~ u 1342 




^yNP 
°1342 


^yNP 
— W231 ) 


^yNP 
u 1423 


,-yNP 
— Ly 2341 




NP 
1342 


^ ^1423 


= 0, 









r<P r»P _ r'NP ^np n 

°1234 ~~ u 2341 — Ly 1234 ~~ u 1423 i u 

/-<NP _ ^NP 
u 1234 — °3421 1 

r NP , r 
^1234 + °: 

(5.18) 

hence only five color factors are linearly independent. The amplitude ( |5.15 ) may be rewritten, for 
example, into the non-symmetric form, 

= g 6 C2341 (^G1234 + ^G234l) + ^1342 (^G1342 + ^4c7342l) + ^1423 (-^G1423 + ^G423l) 

, r NP (A NP , 4 NP 4 NP 4 NP 4 P 94 P 4 P \ fc in) 

' °1423V /1 G1423 "+" ^62341 ~~ ^G1234 ~~ /i G3421 ~~ /i G3421 ~~ Z/1 G1234 ~~ yl G423lJ \ ■>■!■■> I 

, r NP M NP I a NP _ 4 NP _ 4 NP _o 4 P _ 4 P _j_ 4 P 
+ ^1342*. G1342 ' G4231 ^01234 ^G3421 Z/1 G3421 /1 G1234 + ^04231/ 

where each color factor is linearly independent from the others for SU(N C ) with a generic value of 
N c . An advantage of a linearly independent form like eq. ( 5.1SQ is that the coefficient of each color 



factor must satisfy supersymmetry identities by itself. The particular form ( pM9 ) is convenient for 
analyzing the s-channel three-particle cuts shown in fig. |j. 
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Figure 6: Parent diagrams for the fermion loop F contributions. 




Figure 7: Parent diagrams for the scalar loop S contributions. 



We mention that although the original construction [43] of the color organization ( |5.15| ) made 
use of special properties of N = 4 supersymmetric amplitudes, the decomposition actually holds for 
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(Pi 



(Pa 



(P 3 ) 






(NP 2 ) 




^ 3 



(NP 3 ) 



Figure 8: Parent graphs for the mixed M fermion scalar loop contributions. Diagrams Pi, P2, P3, NPi, 
NP2 and NP3 contain Yukawa interactions and can contribute to either the pure gauge coupling case, or the 
superpotential case proportional to |£| 2 ; whereas diagram P4 contributes only to the pure gauge coupling case. 
In the pure gauge coupling case, the solid lines with arrows represent matter fermions and the solid lines with 
no arrows gluinos. In each case there are additional contributions obtained from swapping the scalar lines 
with the matter fermion lines, but this swap does not alter the color factor. 



any purely adjoint two-loop amplitude, as long as all vertices are built out of f abc factors. In such 
cases, one can make use of the Jacobi identity to rearrange the color factor of any Feynman diagram 
into a combination of those appearing in eq. ( 5.15j ) |ft7|] . 



5.3 Fundamental representation two-loop color decomposition 

The color decomposition for four-gluon two-loop amplitudes containing fundamental representation 



particles in the loops is similar to the purely adjoint case (5.13), 



^fund 



9 6 [("^1234 ^X1234 + ^3421 ^X342l) + ^(234) 



(5.20) 



where each term in the sum corresponds to a parent diagram in figs. pHj-tj- I n this case, it is useful 
to include separate S and F contributions with a closed gluino loop, and either a closed matter 
fermion or a closed scalar loop, as depicted in figs. || and [ji]. Again the label X £ {S, F, M, S, F} 
specifies the figure containing the parent diagram. The values of the F®' fundamental representation 
color coefficients are read off from these figures by dressing the parent diagrams with standard 
Feynman rule color factors, as in sect. [iT^ but taking into account the different representations here. 
For example, for the fermion loop F amplitudes, the color factor associated with each primitive 
amplitude is obtained by assigning an f abc to each three-gluon vertex and a (T a )j to each quark- 
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anti-quark-gluon vertex, and summing over the two directions of the fermion arrows 

The color factors obtained by dressing the Pi and P2 diagrams in figs. 
{T ay_ 



factors are not equal to the color factors -^34 or ^234 obtained from figs. |6| 



and [To] with f abc and 
they differ 



in the Tr(T ai T a2 )Tr(T a3 T a4 ) term. Nevertheless, after the permutation sum in eq. (|oT20l) , this 
discrepancy cancels, in the same way that the discrepancy between and C^ 34 cancelled in the 

purely adjoint representation case. This feature allows us to maintain a uniform set of color factors 
for all contributions with fundamental representation matter; i.e., in eq. 
include an X label on the color factors. 



72w there is no need to 
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Figure 9: The F contributions with a closed fundamental representation matter fermion loop and a closed 
gluino loop. The lines with the arrows represent matter fermions. 




Figure 10: The S contributions with a closed fundamental representation matter scalar loop and a closed 
gluino loop. 



For the fundamental representation contributions to the four-gluon amplitudes, the explicit values 
of the color factors are 

F l2U = F xiu = N c [ Tr (T ai T a2 T a ' A T ai ) + Tr(T ai T as T a2 T ai )} + 2Tr(T ai T a2 ) Tr(T a3 T a4 ) , 

F vm = F i2U = — - [ Tr (T ai T a2 T a3 T a4 ) + Tr(T a4 T a3 T a2 T ai )] + 2Tr(T ai T a2 ) Tr(T a3 T a4 ) , 

9 

jP^jj^ = r J 1 r ^J na i J ia 2 2~' a 3 ipa,4^ _|_ r^,^-pa4rpasrpa2rpai^ Tr(T ai T a2 ) r Q'(T a3 T a4 ) 

^1234 = F i2U = 2 Tr (T ai T a2 ) Tr(T a3 T a4 ) , 

pNPj = _ 2 [Tr(T ai T a3 ) Tr (T a2 T a4 ) + Tr(T ai T a4 ) Tr(T a2 T a3 )] . 

(5.21) 

We have dropped terms containing Tr(T a ») from eq. (5.21), since they vanish for SU(N C ) gluons. 
Such terms would have to be restored in order to describe twodoop amplitudes with a mixture of 
external gluons and photons, as in ref. p0| . 

As in the case of the adjoint representation color decomposition ( |5.14| ), the decomposition ( p. 20 ) 
is not unique; there are linear relations between the various color factors. In particular, in the case 
of four external gluons, after eliminating redundant color factors, using the explicit representation 
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of the color factors in terms of color traces Q5.21 ), one finds a total of six independent color factors. 
Also note that the separation into primitive amplitudes is in many cases artificial. For example, 
since F^^ = ^1234 > the distinction between ^1234 an d ^1234 * s mo °t, when all four external legs are 
gluons. 



5.4 Assembly of complete amplitudes 



We can specify the amplitudes in various theories in terms of the primitive amplitudes introduced 
above, and presented in sect. |6| for the ++++ helicity configuration. For example, the four-gluon 
amplitude in QCD with nt massless quark flavors is, 



-4qcd 



a adjoint 



+ n/ ^ nd(1) 



(5.22) 



The pure glue contribution j^i omt [ s given in terms of primitive amplitudes via eq. ( |5.15|) , where 
the sum in that equation runs over the two diagrams in fig. p| Similarly, the single quark loop 



contribution .A^ nd ^ is given via eq. ( 5.2CQ , where the sum runs over all diagrams in fig. |6| except 



fund (2) 



(P4), which gives the double quark loop contribution A 

Using the supersymmetric Lagrangians in sect. ||| as a guide, we may combine the primitive 
amplitudes given in sect. |6| into supersymmetric combinations, which must satisfy the SWI (|2.3| ). 
The amplitude for N = 1 super-QCD with rif matter multiplets in the fundamental representation 
is, 



A 



N=l 



A 



adjoint 
G 



+ A 



adjoint 



_|_ n ^ ^ w 4 f ^ nd ( 1 ) _|_ v 4 f ^ nd ^ 1 ^ + y^ndCl) _j_ ^fund(l) ^fun<li I i 



,fund(2) 



(5.23) 



which is expressed in terms of primitive amplitudes via eqs. ( {Tig) , (|l|) , and (f30|) . The sum over 
parent diagrams of the form A1234 for each contribution is 



a adjoint 


Fig 


|,{A} = 


{P,NP}, 




a adjoint 
J\p 


Fig 


§,{A} = 


{Pi,P 2 ,P 3 , 


P 4 ,NP 1 ,NP 2 ,NP 3 }, 


<fund(l) 

Ap 


: Fig 


•§,{A} = 


{Pi,P 2 ,P 3 


NPi,NP 3j NP 3 }, 


/ ifund(l) 


: Fig 


•0,{A} = 


{Pi,P 2 ,P 3 


P 5 ,NP 1 ,NP 2 ,NP 3 }, 


.fund(l) 

A M 


: Fig 


•!,{A} = 


{Pi,P 2 ,P 3 


NPi,NP 2! NP 3 }, 


« fund(l) 
F 


: Fig 


•i,{A} = 


{Pi,P 2 }, 




/ ifund(l) 


: Fig 


•0,{A} = 


= {Pi,P 2 }, 




,fund(2) 
J\p 


: Fig 


•|,{A} = 


{P 4 }, 




.fund (2) 
A. s 


: Fig 


•0,{A} = 


{P 4 }, 




,fund(2) 
% 


: Fig 


• |,{A} = 


{P 4 }. 





(5.24) 



Note that A f g nd ^ does not receive a contribution from P5 because of a cancellation between fields 
Ai and Ai circulating in one of the loops. 
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For the case of a matter multiplet in the adjoint representation and the superpotential W = 
|<7^Tr$ 3 , the contributions depending on £ are 



/(Yukawa lt|2/ /(Yukawa , /|Yukawa\ ft- nir\ 

Av=i = Isl (As +A M ), (5.25) 
where the contributing parent diagrams are, 

J$*™: Fig. , {A} = {P 5 } , 

^Yukawa. Fig . | , = {p i; p 2> p 3; NPl) N p 2 , N p 3 } . 



In eq. ( |5.25 ), the implicit color factors appearing are generated by the "dressing rules" given in 



sect. |5_2|. 

One can also consider theories with a higher degree of supersymmetry. For example, a pure 
N = 2 super- Yang-Mills theory contains an N = 2 vector supermultiplet composed of a vector, two 
gluinos and two scalar states. This content can be viewed as an N = 1 theory with one adjoint 
matter multiplet. Modifying eq. ( |5.23D for rif = 1, the pure N = 2 amplitude is 

^adjoint _j_ ^adjoint _|_ ^adjoint (1) ^ ^adjoint(l) ^ ^adjoint(l) 

_|_ ^adjoint(l) ^adjoint(l) ^adjoint(2) ^adjoint(2) ^adjoint(2) 
a adjoint . q yiadjoint(l) <adjoint(l) «adjoint(l) . » <adjoint(2) <adjoint(2) „ ,.adjoint(2) 

(5.27) 

In the second step we have used trivial identities for the adjoint representation, 

a adjoint ..adjoint(l) jadjoint(2) 

J\p — J\p T J\p , 

^adjoint(l) _ ^ ^adjoint(2) ^ 2g~j 

/ iadjoint(l) / iadjoint(2) 

As — .A M 

The contributing parent diagrams for ^ oint « 

are the same as for ^ nd(i) in eq. ( |5.24| ), except that 
the associated color factors are the ones for the adjoint representation, i.e., C^za- (For the N = 2 
case, the arrows on the fermion lines in the figures are not important, although in fig. |8| they can be 
used to distinguish between the two species of gluinos.) 



6 Identical-helicity two-loop amplitudes 

In this section we review the results previously obtained for the pure gluon and scalar loop contri- 
butions to the +- 1- helicity amplitude [15]. Then we present our results for the QCD fermion 
loop contributions, followed by the contributions involving Yukawa couplings, which contribute to 
some of the supersymmetry Ward identities. The amplitudes presented in this section have not been 
renormalized; in sect. || we discuss their renormalization. The overall normalization of the primitive 
amplitudes below are such that for the gauge case they can be inserted directly into eqs. ( |5,13| ) and 
( 5.20[ ) without any additional combinatoric factors. c 
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6.1 Pure glue primitive amplitudes 



Using eq. (5.15), the pure gluon two-loop amplitudes are conveniently expressed in terms of planar 
and non-planar primitive amplitudes whose explicit values are |i~5[ ]: 



^Gi234 = P\ S12TX (D s - 2)(A£ X q + Ap X 2 +q + X q X 2 p+q ) + 16f (A„ • A 



y q ) - A p A q 



(si 2 ,S 23 ) 



+ A(D s -2)I t 



bow-tic 



l(X 2 p + X 2 q )(X p -X q )](s 12 ) 



(6.1) 



+ 



7-bow-tie 
' x 4 



^l^q ((P+ qf + S 12) (S12,S 23 , 



\ 2 X 2 X 2 
yq) - A p A q 



(S12,S 2 3), (6.2) 



^G1234 - PS12I4 P (D s - 2)(Ap A 2 . + Ap Ap +(? + A 2 . Ap +(? ) + 16 ((X p • A g 

where the parent diagrams corresponding to these two primitive amplitudes are displayed in fig. ||, 
and p is defined in eq. ( 5.11 ). The other primitive amplitudes appearing in eq. ( 5.1 5| ) are just 
relabelings of these two basic ones. 



The two-loop momentum integrals appearing in eqs. (|6.lD and (6.2) are defined as follows. The 
planar double box integral 

d D p d D q V(Xi,P,q,ki) (6-3) 

(2tt) d (2tt) d p 2 q 2 (p + q) 2 {p - h) 2 (p - h - k 2 ) 2 (q - k 4 ) 2 (q - k 3 - k 4 ) 2 

is displayed in fig. |Tl|(a). The numerator factor V(Xi,p, q, fcj) is a polynomial in the momenta. The 
vectors X p , X q represent the (— 2e)-dimensional components of the loop momenta p and q. We also 
define X p = X p ■ X p > 0, X q = X q ■ X g , and X p+q = (X p + X q ) 2 = X p + X q + 2A P • X q . The bow-tie 
integral z\ ov, ' tie shown in fig. 11 



c) is defined by 



X 



t° w - tie [P(Xi,p,q,h)}(s 12 ) 



d D p d D q 



V(Xi,p,q,ki 



(6.4) 



(2tt) d (2tt) d p 2 q 2 (p - h) 2 (p - k x - k 2 ) 2 (q - k A ) 2 (q - k 3 - k 4 ) 2 
The non-planar double box integral, depicted in fig. |Tl"|(b), is given by 
Xf p [V(Xi,p, q, ki)](s 12 , s 23 ) 



dPp dP q 



V(Xi,p,q,ki) 



(6.5) 



{2tt) d (2tt) d p 2 q 2 [p + q) 2 [p - k^ 2 (q - k 2 ) 2 [p + q + k 3 ) 2 [p + q + k 3 + k 4 ) 2 ' 

The explicit values of the integrals, as a Laurent series in e through 0(e°), and expressed in terms 
of polylogarithms |4S], may be found in appendix A of ref. [15|. We have checked that these values 
agree with results obtained using general integration methods |23|, 24, 25, 26, p?]]. The integral 



rbow-tic 



l([Xl + X 2 q )(X p -X q )](s 12 ), 



(6.6) 



in eq. ( |6.1[ ) vanishes identically, due to antisymmetry of its integrand in X p — > — X p . However, in 
constructing a consistent set of cut integrands it is useful to keep it around, as well as related 
integrals found in the next subsections. 
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Figure 11: Integral topologies appearing in the two- loop identical- helicity amplitudes, and the loop mo- 
mentum routings: (a) the planar double box integral; (b) the non-planar double box integral; (c) the bow-tie 
integral. The arrows here denote the direction of momentum flow. 



6.2 Scalar loop primitive amplitudes 

The planar scalar loop amplitudes have already been presented in ref. [15], where they were used as 
a guide for the construction of the pure glue amplitudes discussed above. Here we reorganize the 
results a bit to be compatible with eq. ( |5,20| ) for the full amplitudes for the case of fundamental 
representation scalars. We also present the non-planar contributions. 

The scalar loop primitive amplitudes corresponding to the planar parent graphs in fig. ^ are, 



^51234 - 2 pj s 12 2f [Ap Ap +(? ](S12,S23) 



«12 



+ 4J bow-tie [A 2 (Ap . A(?)](si2) | 
^S1234 = 2 p|si2 2"I[Ag A2 +g ](si2,S23) 



+ ^±^lt^[Xl X 2 ((p + qf + , 12 )]( S12 , S23 ) (6 - 7) 
+ 4J bow-tic [A 2 (Ap . A(?)](si2; 

4i234 = 2pl [Sl2 l![X 2 p X 2 q ](s 12 ,s 23 ) + ^™- tic [X 2 p X 2 q ]( Sl2 )} 

A p ((P + + h S ^)\ ( S 12, S23) , 
^1234 =^4 b ° W - tle ^A2]( Sl2 ). 



Ap 

•7-bow-tie 

/1 S1234 — x 4 
•S12 

bow-tie r \ 2 \2i 



Similarly, the scalar loop primitive amplitudes corresponding to the non-planar parent graphs in 
fig. are, 

A S1234 = ^PSl 2 l4 [\ Ap +q \[Sl 2 , S 23 ) , 

^5i P 234 = 2p Sl2 I^[X 2 q A^ +g ](si2, S 23 ) , (6.8) 
A^ 2 % 4 = 2ps 12 ir[X 2 X](s 12 ,s 23 ). 

The integrals appearing in eqs. ( |6.7| ) and (^) appear in the pure gluon case as well; see appendix A 
of ref. [15] for their explicit values through O(e ). We remind the reader that the number of scalar 
states propagating in the above primitive amplitudes is 2 x (N 2 — 1) in the adjoint representation case 
(when they are dressed with color using eq. ( 5.13| )), and 4xiV c states for the N c + N c fundamental 
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representation case (when they are dressed with color using eq. ( |5,20| )). In either case, this number 
matches the number of states of the fermion superpartners in the fermion loop primitive amplitudes 



presented in sect. p. 5 . 

A convenient way of organizing the numerators of the planar and non-planar double box inte- 
grals is in terms of the labels of particles circulating in the loops. Because the external momenta 
ki,k 2 ,k 3 ,k4 are defined to have vanishing extra dimensional components, as one follows a scalar 
line around only two different extra dimensional momenta appear, which we can label by A Sl and 
X S2 . For example, in fig. 0(Ps) the two distinct extra dimensional scalar momenta are \ p and X q . 
Using these labels, the numerator arguments of the planar and non-planar double box integrals in 
eqs. ( |6.7D and (|6.8|) all have a uniform structure, 



\ 2 \ 2 



(6.9) 



This observation simplifies the bookkeeping when evaluating the three-particle cuts. 



6.3 Fermion loop primitive amplitudes 



The structure of the fermion loop primitive amplitudes is quite similar to that of the scalar loop 
primitive amplitudes. The parent diagrams are described by fig. ||, and the results are, 



A 



F1234 



A P2 

^F1234 



A 



F1234 



^F1234 



1NP1 
1 F1234 

iNP 2 
1 F1234 

,NP 3 
1 F1234 



p\ s i2^4 ~\{Ds - 2)(A p+ g • A P A 2 ) - 2Xp\p +q - 4((A P • X q ) 2 - X 2 X 



2\2^ 



,Sl2,S23j 



- 2 



.D,-2 



^bow-tie r \ 2 \ 2 



S12 



4 l"p"q 



[X 2 X{(p + qy + 51 2 )](S12,S2 3 ) 



l "-- tie [A^ p -AJ( Sl2 ) 



p\ si 2 ll\-\{D s - 2)(X p+q ■ X q X 2 p ) - 2X 2 q X 2 p+q - 4((Ap • A g ) 2 - X 2 p X 2 q )} (s 12 , s 23 ) 



- 2- 



D x 



rbow-tie r \ 2 \ 2 



S12 



(D s + 6)l^- tie [X 2 X p -X q ](s 12 : 



p\ 812 ll\\{D a - 2)(A P • X q X 2 p+q ) - 2X 2 p X 2 q - 4((A P • A 9 ) 2 - X 2 p X 2 )} (s 12 , s 23 



+ (D s -2)-2 )s l2 Z 



r-bow-tie 



bow-tie 



(X p -X q )[2(X z p + X 2 q 



[X p - X q ](s 12 ) + (D S -2)X 



bow-tie r \ 2 \ 2 



5S12] 



+ 



312 



•7-bow-tie 
--t-4 



ApA 2 ((p + g) 2 + \si 2 ) (S12,S 2 3 



(6.10) 



p Sl 2lf P [-^(D s - 2)(X p+q ■ XpX 2 ) - 2X 2 p X 2 p+g ~ 4((Ap • A g ) 2 - A 2 A 2 )]( Sl2 , S23 ) , 
, = ps X2 lf > [-\{D s - 2)(X p+q ■ X q X 2 ) - 2X 2 q Xl +q - 4((Ap • X q ) 2 - A 2 A 2 )](s 12 , s 23 ) , (6.11) 
A 7i2M = psi 2 lT[\{Ds - 2)(X P ■ X q X 2 p+q ) - 2A 2 A 2 - 4((A P • A,) 2 - A 2 A 2 )](s 12 , s 23 ) . 
As was the case for the scalar loop amplitudes, the normalization of the fermion loop amplitudes 
(when dressed with color using eqs. ( 5.13| ) and (5.20)) corresponds to 2 x (A" 2 — 1) states for the 
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adjoint representation case and 4x N c states for the N c + N c fundamental representation case. Most 
of the integrals in eqs. ( 6.10| ) and ( |6.11 ) are again the same ones as for the pure glue case [15]. There 



are, however, a few new integrals. In ref. |19[ these amplitudes, together with the other helicity 
configurations, are given directly in terms of polylogarithms, after subtracting off certain universal 
pole terms in e pC| ]. 

Just as we did for the scalar loop amplitudes, we can write the numerator arguments of the loop 
integrals in terms of the two distinct extra-dimensional momenta flowing in the fermion propagators, 
which we call A/ x and A/ 2 , and the one in the internal gluon propagator, labeled as X g . In this case, 
the numerators appearing in the planar and non-planar double-box integrals all have the form, 

-\{D S - 2)(A /1 • A /2 A 2 ) - 2X\Xl - 4((A /x • \ h f - X%X%) , (6.12) 

using the fact that 

(A p • X q ) 2 - X 2 X 2 = (Ap+q • A p ) 2 - Ap +(J Ap = (X p+q ■ X q ) 2 - X 2 p+q X 2 q . (6.13) 
6.4 Mixed scalar and fermion amplitudes 

In order to fully study the supersymmetry Ward identities, we need also the "mixed" contributions 
where both fermions and scalars appear in the loops. The parent diagrams associated with these 
contributions are shown in figure |8[ In these diagrams, the solid line represents fermions, while the 
dashed line represents scalars. These diagrams, and primitive amplitudes, describe the pure gauge 
coupling contributions A^ A , for which one of the fermions is a gluino, via eq. ( |5.24 ). They simulta- 



neously describe the superpotential (^-dependent) contributions, for which both of the fermions are 
matter fermions, via eq. ( 5.26) ). 



Following the notation of the previous sections we find that the explicit forms of the planar mixed 
primitive amplitudes are 

^M1234 = p{ s 121a[K ■ \^p+q ~ \ ' \+q^p\ ( s 12, S 23 ) + 512X4°™ tlG [A p ■ A g ](si 2 )} , 

4)1234 = p{si2l![X p ■ X q X 2 p+q - X p ■ X p+q X 2 q ](s l2 , s 23 ) + s 12 l^- tic [X p ■ X q ](s 12 )} , (6.14) 

^M1234 = -p{ s 122f[Ap • X p+q X 2 + X q ■ X p+q X p ]{si 2 , S 23 ) + 4l\° W tl6 [X p X 2 ] (s i2 ) } , 

corresponding to the first three diagrams in fig. ||. For the case where the fermion and scalar loops 
are separated from each other, which does not contribute in the superpotential case, we have 

^1234 = ^{xr^IAp • X q {X 2 p + X 2 q )]{8 l2 ) + i-Tb°w-tic [ A 2 A 2 ((p + g) 2 + 1^)] ^ ^ 

(6.15) 

The mixed contributions corresponding to the three non-planar parent diagrams in fig. || are 

^M1234 = P s l2^f P [X p ■ XqX^ +q - X q ■ X p+ qX p ](si 2 , S 23 ) , 

^M1234 = P s i2^i P [X p ■ XqXl +q - X p ■ X p+q X 2 ](si 2 , S 23 ) , (6.16) 

^M1234 = -P s 12^f P [X p • X p+q X 2 + X q ■ X p+q X 2 )(si 2 , S 23 ) . 
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As in the previous cases, in the double box integrals we can label the numerator momentum factors 
in terms of which types of particles carry the loop momenta, i.e., 

-X h -X h X 2 s . (6.17) 

In this case there are two contributions to each primitive amplitude above, since interchanging matter 
fermions with scalars in fig. |8| does not alter the color factor. 



6.5 Amplitudes with both a gluino and matter loop 

Up to color factors, the parent diagrams with both a gluino and matter loop depicted in figs. || and 
|l0| are identical to the parent diagrams in fig. ^(P^) and fig. |8|(P4). Thus, the associated primitive 
amplitudes are also simply related, 

A Pi _ 4P2 _ 4P4 

^F1234 ~ A F1234 ~ A F1234 . 

4x234 = -2p{x 4 bOW " tiC [A p • X q \ 2 q ](s 12 ) + Ajbow-tic [ A 2 A 2 ((p + q) 2 + 1^)] ^ ^ ^ 
4l234 = "2p{x 4 bOW - tiC [A p • A 9 A^]( Sl2 ) + Ajbow-tie [ A 2 A 2 ((p + ^2 + 1^ 



V S12,S 23/ 



Hence 

1P2 



A §1234 + ^51234 ~ A A/1234 > ( 6 - 19 ) 

icorporates the 

in diagram (P4) of fig. pi while the diagrams in fig. [Hi separate the two cases. 



accounting for the fact that ^^1234 incorporates the cases were the scalar loop is on the left or right 



7 Two-loop supersymmetry identities 

We now verify that the supersymmetry Ward identity fl2,3| ) does in fact hold for the amplitudes 
presented in sect. || when using the FDH scheme. 



7.1 N = 1 Identities 

The two- loop N = 1 SWI are obtained by applying eq. (2J3) to eqs. (|5.23j ) and ([5.25 ). For conve- 
nience, we subdivide the SWI according to the independent parameters appearing in the amplitudes, 
rif and £, 

^adjoint + ^adjoint = Q ^ 
^fund(l) _j_ ^fund(l) ^fund(l) ^fund(l) ^fund(l) _ q ^ ^ 

A^ nd{2) + ^ und{2) + ^ nd(2) = , (7.3) 

^Yukawa + ^Yukawa = q ^ ^ 

corresponding to the four sets of diagrams represented in fig. ^. These identities must hold if 
supersymmetry is preserved by a regularization scheme. 
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In general, we can further subdivide the identities by considering separately the coefficient of each 
linearly independent color factor. It turns out, however, that for the identical-helicity amplitude a 
slightly stronger subdivision is possible directly in terms of the primitive amplitudes presented in 
sect. ||. In other words, we shall show that the identities hold for each coefficient of the different 
color factors appearing in eqs. (5.14) and ( 5.20| ), even though they are not all linearly independent. 

The following two integrals, which will appear on the right-hand-side of several of the identities, 
have a simple representation in terms of double-box integrals evaluated in 6 — 2e dimensions, 



K p = ps 12 l P {\ p -\ q f 



\ 2 \ 2 



-rjNP _ n a tNP 



(S12,S23) : 
(A p • \ q f - ApA^ (S12,S 23 ) 



; (l + 2e) (47r) 2 ps 12 ll 
J(l + 2e) (4ir) 2 ps 12 l 



-2e, 



lSl2,S 2 3j 



NP,D=6-2e 



(S12,S 2 3) 



(7.5) 



where 



T P,D= 


(Sl2,S 2 3j 

r d 6 ~ 2e P 


d 6 ~ 2e q 






1 




~J (27r) 6 - 2e 


(27r) 6 - 2e 


p 2 q 2 (p + t 


i) 2 (p-h 


)2 ( p - kl - k 2 ) 2 (q - k 4 ) 2 (q-k 3 - k 4 ) 2 ' 




(•Sl2,S23j 

r d 6 - 2e p 


d 6 - 2t q 






1 




~ J (2vr)6-26 


(2^) 6 " 2e 


p 2 q 2 (p + q 


) 2 (p - h 


) 2 (q ~ k 2 ) 2 (p + q + h) 2 (p + q + k 3 + k 4 ) 2 



(7.6) 

The planar and non-planar double box integrals have neither infrared nor ultraviolet divergences in 
six dimensions. Therefore, due to the explicit e in front of the integrals in eqs. (ffl]), 1Z P and 7£ NP 
are both of 0(e) as e — ► 0. 

We begin our inspection of the supersymmetry identities with eq. (7.1). This identity holds in 
the FDH scheme (D s = 4), because the combinations 



A p 

/1 G1234 



4 NP 
yl G1234 



8=1 

3 



p* 

F1234 



67^ 



NPi 
F1234 



6K 



NP 



0(e) 



0(e) 



(7.7) 
(7.8) 



i=i 



vanish as e — > 0. The right-hand sides in these equations are obtained by inserting the explicit values 
of the primitive amplitudes given in sect. [6|. We dropped the integrals that vanish identically, i.e. 
the bow-tie integrals that are odd in A p or X q . The same equations ( f7v7| ) and (7.8) obviously also 
hold after performing any permutation of the external legs. 



Similarly, the identity (|7.2|) holds because 

2 

I Pi 



E( A Pi 4- /T* 4- A 
V ^F1234 ^ ^ 1 S1234 ^ ^M1234 



i=l 



A P3 -4- V 4 Pl +4 
^\F1234 ^ ^51234 ^ Ji 

ie{3,5} 

2 



p 3 

M1234 



-2TV 



|NP, 
■\F1234 



i=l 



^ /1 S1234 ^ ^M1234 J 



-All 



NP 



0(e) 



0(e) 



0(e) 



(7.9) 
(7.10) 

(7.11) 
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^F1234 + Asi234 + ^M1234 = ~ 2 ^ NP = C(e) , (7-12) 
A Il234 + A Jl234 = ^Il234 + ^1234 = ' ( 7 - 13 ) 

and the identity ([7^) holds because 

^F1234 ~~ Asi234 = ^ ' (7-14) 
2 ^F1234 + ^M1234 = ■ (7-15) 

The latter three equations hold to all orders in e. 

Finally, eq. ( [7.4|) , the identity for the amplitudes depending on the Yukawa coupling £, holds 
because 

I E 4li234 + 2^ 2 34 = 2TZ P = 0(e) , (7.16) 
z i=i 

3 



^E< P i234 = 2^ NP = 0( e ), (7.17) 



where the factors of 1/2 and 2 on the left-hand side of the equations are relative combinatoric factors 
compared to the gauge case. 

The same relations hold for any permutation of the external legs. We conclude that for D s = 4 



(FDH scheme) the identical-helicity (++++) primitive amplitudes satisfy the SWI (2^) as e — > 0, 
at least through two loops. In the course of computing the two-loop amplitude for the helicity 
configuration — h++ in QCD flT9| , we have also verified that the pure super- Yang-Mills identity (7.1 



is satisfied for — h++ up to 0(e) corrections, again provided that D s = 4. For either of these two 
helicity amplitudes, in the HV scheme where D s = 4 — 2e the identities ( \7.T\ ) and ( |7.8| ) with gluon 
loops are not satisfied, even at C(e _1 ). 

Even in the FDH scheme, the above two-loop N = 1 identities are generally satisfied only for 
e — > 0, because the quantities 1Z P and 7£ NP are known at order e, and they are nonvanishing at 
this order. If one thinks about constructing three-loop amplitudes via their unitarity cuts, which 
include the product of a two-loop amplitude with a tree amplitude, one might be concerned that the 
0(e) breaking of the identities at two loops could lead to a non- vanishing breaking of the SWI at 
three loops even as e — > 0. Of course there will also be contributions to the product of two-loop and 
tree amplitudes from intermediate momenta in (— 2e) dimensions, which may well cancel the above 
contributions. Clearly this issue warrants further investigation. 



7.2 Exactness of iV = 2 Identities 

An 0(e) breaking in N = 1 supersymmetry Ward identities in the FDH scheme is perhaps to be 
expected since the scheme is defined by "dimensional expansion", i.e., D > 4. Supermultiplets 
become larger as D increases; hence for D > 4 an N = 1 supersymmetric theory in D = 4 does not 
provide enough states to form a supermultiplet for the intermediate unobserved states with momenta 
in the extra (— 2e) dimensions. Thus, one might expect some kind of violation of supersymmetry 
for finite e. On the other hand, the N = 2 vector multiplet, consisting of a gluon, two gluinos 
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and a complex scalar, may be viewed as coming from the compactification of a higher-dimensional 
supersymmetric theory, for example N = 1 supersymmetric gauge theory in D = 5. So, under a 
dimensional expansion of the N = 2 theory away from D = 4 there are sufficiently many states 
to form a supermultiplet, and we might expect the SWI to hold to all orders in e in the N = 2 
case. The supersymmetry identity that should be satisfied in pure N = 2 gauge theory is given 
in eq. ( |5.27| ). This identity does hold to all orders in e, because of the following relations between 
primitive amplitudes: 

5 3 



^G1234 + E Asi234 + E( 2 ^F1234 + A 



M1234) + 4^F1234 + 2 ^Afl234 



0. 



i=l 



i=l 



which is the sum of eqs. ( |7.7j ), ( |7.9| ), and ( |7. 1Q| ) ; and 



(7.18) 



A 



NP 

G1234 



51234 ^F1234 ' ^M1234^ 



0. 



(7.19) 



i=l 



which is the sum of eqs. (|7j), ( ffHD , and ( |7lp . 

Now consider adding nj hypermultiplets in the N c + N c representation to the N = 2 gauge 
theory. The N = 2 hypermultiplets have the same field content as the pairs of N = 1 chiral matter 
fields, Qi + Qi. The nj terms in the SWI for D s = 4 clearly hold to all orders in e, by virtue of 
eqs. ( 7.14 ) and ( 7.15 ). The same statement is true for the nj terms, however here the cancellation is 
more subtle. One can split the N = 2 gauge multiplet into an N = 1 gauge multiplet plus an N = 1 
adjoint matter multiplet <£. There is a set of nj contributions from coupling the hypermultiplet fields 
to the N = 1 gauge fields, which is given by precisely the combinations (|7.9| )-( 7Tl2| ) encountered 
above. A second set of nj contributions comes from coupling the hypermultiplet fields to <£. This 
coupling is through an N = 1 superpotential term, W oc J2i Qi^Qi- The strength of W is dictated 
by the gauge coupling, because the Yukawa interactions with the N = 1 gaugino A in eq. (|2.5|) have 
to be duplicated by couplings to the N = 1 matter adjoint fermion, which is the second gaugino 
of N = 2 supersymmetry. These Yukawa interactions just double the mixed fermion-scalar terms, 
^Mi234 and A^ 234 . However, there are additional mixed terms due to the Yukawa vertices from 
W which contain the adjoint scalar particle. These terms are dressed a bit differently with color 



than the mixed contributions considered in sect. 5.3. However, they still can be written in terms of 
4-M1234 and A^ l 234 . Finally, there are new P5 terms from scalar four-point interactions generated 
by W . Adding all the terms together, we find that the SWI hold to all orders in e by virtue of 



2 



i=l 



^F1234 ^ ^-51234 ^ ^A/1234 



+ 



3 

E ^ A/1234 + 4^S1234 I 
i=l / 



4^ 3 4- 
^1234 ^ 



E 

i€{3,5} 



^51234 



4- 4 r3 
^ ^M1234 I 



2 



8=1 



NPi 



, 4NP, A 
F1234 ^ ^51234 ^ ^A/1234 



.NP, 



+ 



1 

+ 2 



3 



3 

^Af 1234 ^ ^51234 



0. 



4NP3 , AVir 3 

F1234 ' 1 



|NP 3 
1 S1234 



4NP3 
^A/1234 



+ 



/1 M1234 



3 \ 

V 4 NPi 1 
/ , ^M1234 I 



(7.20) 

(7.21) 
(7.22) 
(7.23) 



which follow from adding appropriate pairs of eqs. (|7.9|) ( |77T7 ) . 
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8 Renormalization 



The identical-helicity amplitudes presented in sect. |6| have not been renormalized. Since the corre- 
sponding tree-level amplitudes vanish for this helicity configuration, the renormalization procedure 
is similar to that for a typical one-loop amplitude. In particular, the modified minimal subtraction 
counterterm to be subtracted from A\ oop (l + , 2 + , 3 + , 4 + ) is 



4g\c r -Al- loop (l + ,2 + ,3 + ,4 + ), (8.1) 



e 

where 



r(i + e)r 2 (i- e) 

T(l -2e)(47r) 

and &o is the one-loop /3-function coefficient. For example, in QCD with N c colors and rif quark 
flavors, 

„«™ = H^3£, (8.3) 
while for pure N = 1 super- Yang-Mills theory, 

tf-^ftfc (8-4) 

More generally, the counterterms depend on the scheme choice through the scheme dependence of 
the one-loop amplitude. 

8.1 Coupling constant shift 

When comparing physical results obtained in the FDH scheme to ones obtained in the 't Hooft- 
Veltman (HV) scheme, one must account for the shift in the coupling. In general, the coupling 
constant of the theory is scheme dependent. The one-loop relation between the QCD coupling 
constant in dimensional regularization and that in dimensional reduction has been computed pre- 
viously g, 0> 0- Here 

we shall give the relation to two loops. We also allow for an arbitrary 



value of the variable controlling the number of gluon states in loops, D s = 4 — 2e<5^, so that one 
can consider schemes other than the FDH/DR scheme (5r = 0) or the HV scheme (5r = 1) if one 
desires. 

A convenient method for calculating the scheme dependence of the coupling constant is to eval- 



uate the shift in the background field gauge method 51, [3^, 52, 23 by computing the vacuum 



polarization, following the one-loop discussion of ref. [11]. The background field gauge possesses 



— 1/2 

a Ward identity Z g = Z A which allows one to obtain the coupling constant renormalization Z g 
from the wave-function renormalization Za- Our two- loop computation follows closely the ones in 



refs. [52, 53|, except that we have to keep track of the finite terms, not just the poles in e, when the 
terms are scheme dependent. 

The coupling constants of the 't Hooft-Veltman and CDR schemes are identical because the 
vacuum polarizations are computed using exactly the same rules for internal states. The HV and 
CDR regularization schemes, together with the widely used modified minimal subtraction (MS) 
renormalization scheme, define the coupling constant a^jg- The modified minimally subtracted 
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versions of the FDH and DR schemes, which we call FDH and DR, define the coupling constant 
a FDH = q dr- F° r an y coupling a, let the reduced coupling be a = a/(2ir). Then the one-loop 



coupling constant relation is M, 31, 32, [TT| 



(, . Ca 

a DR ~~ a FDH ~~ a MS V 6 Qms ' 

where Ca is the quadratic Casimir in the adjoint representation, Ca = N c for gauge group SU(N C ). 
Our goal here is to compute the two-loop corrections to this relation. 

A priori, there is a difference between FDH and DR, in that DR divides a gluon Lorentz index 
into a D = 4 — 2e index and a residual 2e index. The latter components are referred to as e-scalars. 
Imposing D-dimensional gauge invariance still permits the e-scalars, and their associated couplings, 
to renormalize differently from the D-dimensional gauge fields. For example, when fermions are 
present, their coupling g to the gauge fields gives rise to a Yukawa coupling g of the fermions to the 
e-scalars. In a non-supersymmetric theory, when one also requires the e-scalar Green functions to 
be finite, one finds that g and g have different (3 functions even at one loop |]]. 

A nice feature of the background field method is that none of the quantum fields have to be 
renormalized [p2| . Quantum fields do not appear on the external (background) legs, so Z factors 
cancel between propagators and vertices. Therefore, one does not have to worry in the present 
calculation about whether the e-scalars that are present in DR should be renormalized differently 
from the ci-dimensional gauge fields. Thus, the counterterm graphs are identical between DR and 
FDH. The non-counterterm graphs also turn out to be identical; the D and 2e ranges of the indices 



in DR combine to give exactly the four-dimensional range used in FDH. Thus the FDH and DR 
gauge couplings are equivalent, at least in any theory (supersymmetric or not) with only a gauge 
coupling present. 

Following the computational rules for the FDH scheme, we write the Lorentz contraction that 
tracks the number of physical states as 

riZ = D s = 4-2e5 R , (8.5) 



where r\ is the Minkowski metric, 5r = 1 for MS, 5r = for DR or FDH. In ref. p2| , p3| ] the divergent 
terms in the two-loop vacuum polarization were computed in MS. Here we also need the portion of 
the finite terms that depends on 5r. For pure gauge theory, Table 1 of ref. j52| gives the divergent 
parts of the Feynman graphs shown in fig. |l^. Of the graphs a-m, only b, h, k, and the counterterm 
graphs 1 and m, are 5r dependent, although a little inspection is necessary to show that graphs d 
and e are independent of 5r. 

The counterterm contributions 1 and m come from one-loop renormalization of the gauge-fixing 
terms. For our purposes, this renormalization has to be performed including the finite level, whereas 



in obtaining the f3 function in ref. [52] it was sufficient to carry out the renormalization only through 



0(l/e). We find that the gauge- fixing renormalization factor Z a is given by 



1 + 



5 5 R 28 \ / 4 20 \ 1 g 2 



(8.6) 



(4tt) 2 ' 

where the fermion generators have trace Tr(T a T ) = Tp5 ab (for QCD with n/ flavors, Tp = \ n f)- 
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(c) (d) (e) 




Figure 12: Feynman diagrams for the pure glue contributions to the vacuum polarization, following the 
labeling of ref. J52j. Dashed lines represent ghosts, and square dots represent counterterms. 



In the notation of refs. 



the vacuum polarization is 



n^(fc) 



ig A C\8 ab 



(4tt) 4 



Ag^u Bk^k u 



(8.7) 



In the background field method the final sum over diagrams gives A = B, ensuring the transversality 
of the vacuum polarization. Here we give just the B terms, dropping for convenience the combination 
7_B — ln47r + ln(/c 2 //i 2 ), as well as all Jfj-independent finite contributions. For the B terms of eq. (8.7) 
we obtain, 



graph b: 
graph h: 
graph k: 
graphs 1+m: 
a-m total: 



25 

6^V 
_9 
~ 8?V 
27 

8?^ ' 54 

10 2 . 

T" + T°R > 
3e 3 

17 7, 



49 



2c 



(l + -e 



4e 



(S 



R 



4e 



(Sr - 1 



6 



The values of the remaining graphs with no 5r dependence may be obtained from ref. p^] . 

For the fermion-loop contributions shown in fig. Il3], we find that the only finite (^-dependent 
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(a) 



(e) 



(b) 



(d) 



(f) 



Figure 13: Fermion loop contributions to the vacuum polarization, following ref. 
fcrmions, and square dots represent counterterms. 



Solid lines represent 



terms that survive are in the T F C F color structure, and the total in Table 1 from ref. |53J becomes 

total: - - {—T F C A + 2T F Cp) - 2T F C F 8 R . (8.9) 

Here C F is the quadratic Casimir in the fermion representation, C F = (N 2 — 1)/(2N C ) for the 
fundamental representation of SU(N C ). It is rather simple to see that adding scalars to the theory 
does not induce any additional 5r dependence; the corresponding diagrams do not contain any 
contractions. 

The background field Ward identity [52] relates the inverse coupling directly to the vacuum 
polarization. Using this relation, we find that the bare inverse coupling is given in terms of the 
renormalized coupling a$ R = as R /(27r) by 



^bare 



— \ 1 - a Sf 

a &R 



1(11 

2e V 3 



Ca - -T F 



\c A 5 R + ci 
o 



/e^ E k 2 y 
\ Air ) 



-(-ci 

8eV 3 A 



20, 

T 



T F C A - AT F C F ) + (h& A - ^T F C F )5 R + c 2 ( 



4tt 



_ (8.10) 

where c\ and c<i are independent of 5 R . To relate as R to the standard MS coupling, a = aj^g/(2-7r) 
with 5 R = 1, we equate the two versions of the bare inverse coupling, which leads to 



1 + ho A {\ -5 R )a+ (^Ci(l - 5 R ) + ^C|(l - 5 R f - l -T F C F {\ - 5 R ) 



$.111 



In particular, for a gauge theory with a general fermion content, the DR coupling, a = a^^/(2iv), is 
expressed in terms of the MS coupling by setting 5 R = in eq. (8.11), yielding 



a = a 



1 '" - l T F c F y* 



1 + - R c ^ + (th c - 



18~ A 2 



+ 0(a 3 ) 



.12) 



8.2 Three-loop (3 function 

An interesting side-benefit of the above computation is that it allows us to obtain the three-loop (5 



function in the FDH and DR schemes from the known (3 function in the MS scheme [33]. 
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In general, knowledge of the scheme dependence of the coupling allows one to convert the (5 
function between schemes at one higher loop order. For example, the /3-function coefficients in 
four different momentum-subtraction (MOM) schemes were computed at four loops in similar fash- 
ion [55]; see also refs. [56] for three-loop MOM computations. There are similar computations on 
the lattice [57]. 

Consider two different renormalization schemes, and hence two different coupling constants, for 
a single theory, 

(8.13) 



o 



a 



2vr ' 2vr 
and suppose they are related through two loops by 



a = a 
a = a 



1 + d\a + d 2 a + 



l-d 1 a + (2d'( - d 2 )d 2 + . 



U4) 



Let the functions in the two schemes be 

da „, , 

^ = «•) = -» 

pf. ^ 0(H) ee -2 
djj, 



b a 2 + b ia 3 + b 2 a A + ... 



b d 2 + bid 6 + b 2 d" + . . . 



.15) 



One can also calculate the function in the second scheme as 

da da dd s r _ „ , 9 , 

/3(a) = Ht- = x T" = 0( a ) x I 1 + 2d i a + M 2 a 2 +•••]• 
dfj, d\i da 



(8.16) 



Comparing this expression with that in eq. ( B.15| ), after substituting for a using eq. ( |S.14 ), we have 

[b a 2 + bia? + b 2 a 4 + . . .] x [1 + 2d x a + 3d 2 a 2 + . . .} 
= b a 2 [l + d ia + d 2 a 2 + . . .] 2 + 6 ia 3 [l + d x a + . . .] 3 + b 2 a 4 [l + . . .] 4 ' 

Equating coefficients of a gives the standard results that bo = bo and b\ = b% , and also 



3.17) 



d-,b - d\b\ + d 2 b . 



.18) 



This gives the shift in the three-loop /3-function coefficient in terms of the one- and two-loop coupling 
constant shifts d\ and d 2 . 



To compute the three- loop /3-function coefficient in DR scheme, we now use eq. ( $.18 ), in con- 
junction with the well-known three-loop /3 function in MS [p£ 



j nil, 



fro = -^Ca - rl> , 
6 3 

17 5 
b\ = —C A — -CaTf — CfTf 



H = kid 



3.19) 



1415 



Cirri 205 ± 2 rri 1 ^ /~1 rri2 , X X *~i rjil 

A 1 F T^^A^F^F + -CplF + Tno^A 1 F + To^F 1 F 



432 A 216 
Hence the three-loop /3-function coefficient in DR is 

1 



79 
108 



11 
— ( 

18 



79 



17 



b 2 — C\ C 2 aTf C aC fTf H — C%Tf H CaT'w H CfTj? . 

432 A 216 A 72 4 F 108 F 18 F 



(8.20) 
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For QCD, SU(3) with riy flavors of quarks, the result is 

8.21 



3115 5321 373 2 

nt H n f . 

16 144 1 432 ' 



For a generic number of flavors, this value is significantly closer to the MS value than are the MOM 



versions of b 2 reported in refs. [56, although for five flavors the value in the scheme labeled 
MOMggg in ref. |56| happens to be the closest to MS. The result ( |S.20| ) should be obtainable as an 
extension of the work of ref. |54j . In that paper, three-loop f3 functions were calculated for a variety 
of theories, in a way that could be applied to both MS and DR schemes, but only the MS results 
were presented. 

For the special case of pure super- Yang-Mills theory, with a single Majorana gluino, we set 
Tp = Ca/2, Cf = Ca, to obtain 

3„ ,3 



19 r3 r 21 3 



(8.22) 



in agreement with ref. ||6^] . (The MS value, b 2 , was first calculated in ref. [pq].) 
8.3 Scheme dependence of the primitive amplitudes 

Since our expressions for the amplitudes in sect. ^ smoothly interpolate between the FDH (D s = 4) 
and HV (D s = 4 — 2e) schemes we obtain the scheme dependence of the amplitudes by subtracting 



the two forms. First consider the pure glue amplitude. From eqs. ( |6.lD and (6.2), that the difference 
between the unrenormalized pure glue primitive amplitudes in the two schemes is: 

° A GX2U = ^G1234 ~ A G12M = ~ 4 — ~ ~ X A [\\ + 

) tie (8-23) 
= " 2 7 ^ 1U ° n loop ' FDH (l 5 + , 2+ 3+ 4+) + 0(e) , 

where we used eq. (A. 11) of ref. [15] to obtain the first line and eq. ( [5,10| ) to obtain the second. 
Similarly, for the non-planar pure glue primitive amplitudes 

SA^U - <234 DH " «4 V = ^ > ^ , 3+ 4+ ) + <>® > (8-24) 

where we used eq. (A. 14) of ref. J15|] . In these expressions, the coupling associated with the HV 
amplitude is the standard MS one, and the coupling associated with the FDH amplitude is the DR 
one. This shift does not include the coupling shift ( 8.11] ) or the shift in the counterterm, which are 



straightforward to incorporate. For the remaining cases with fermion, scalar, or mixed loops we find 
that the differences are all of 0(e). For example, 

<54i234 = -ep{s l2 li[\l\ 2 q + X p ■ X q X 2 q ] + 2T 4 bow " tic [A^(A P • X q + ^-X 2 p ((p + q) 2 + s l2 ))] j = 0(e) , 
^234 = -eps 12 lf p [^ 2 q + A P • X q X 2 q ] = 0(e) . 

(8.25) 
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The simple structure of the primitive amplitude scheme shift is, of course, not accidental. It may 
be understood in terms of the infrared divergences |50| appearing in the amplitudes. As one shifts 
between the schemes the integrands undergo shifts of 0(e). These shifts can contribute only when 
they are multiplied by e -1 divergences. The complete scheme dependence through 0(e°) therefore 
depends only on the scheme dependence of the universal divergences. The scheme shift for two- 
loop helicity amplitudes for which the tree amplitudes vanish are especially simple, because the 
divergences are rather like those of typical one- loop amplitudes. One may also shift the amplitudes 
of sect. |6|to the CDR scheme following the same strategy used at one loop [11, 30 1. For more general 



helicity amplitudes, the scheme shift is more complicated because one must account for non-trivial 
scheme dependence in the /^(e, fj, 2 ; {p}) function of Catani [50|. These issues are presented in more 
detail in ref. (lj| . 

9 Conclusions 

In this paper we presented detailed rules for the four-dimensional helicity (FDH) regularization 
scheme to higher loops. The scheme is designed to preserve the number of bosonic and fermionic 
states at their four-dimensional values. Computationally, it amounts to a relatively minor modi- 
fication of the 't Hooft-Veltman |2| dimensional regularization scheme. As an initial check of the 
supersymmetric properties at higher loops, we explicitly evaluated a number of two-loop exam- 
ples, illustrating that the scheme preserves the required supersymmetry Ward identities, while the 
't Hooft-Veltman scheme fails (as expected) to do so. The identities checked in this paper were 
relatively simple to evaluate directly, because they involved identical-helicity four-gluon amplitudes 
which have an especially simple analytic structure. Although we did not present it here, we verified 



that for N = 1 super- Yang-Mills theory, the — h++ helicity amplitude [19| also vanishes in the FDH 
scheme, in accordance with supersymmetry. It would be interesting to systematically investigate the 
supersymmetric properties of the remaining helicity amplitudes. In these cases the Ward identities 
relate amplitudes with four external gluons to amplitudes with external fermions as well; the former 
have been computed |l5| but the latter have not. 

By keeping track of a parameter 5r which smoothly interpolates between the 't Hooft-Veltman 
and four-dimensional helicity schemes, one can conveniently verify that supersymmetry Ward Iden- 
tities hold even in non-supersymmetric theories such as QCD — after modifying color factors and 
multiplicities so that states fall into supermultiplets. This approach is of some interest in non- 
trivial QCD calculations, where it can serve as a cross check. Indeed, the fact that the two-loop 
identical-helicity amplitude presented here satisfies the required supersymmetry Ward identities pro- 
vides one additional check on the four-gluon QCD amplitudes of refs. ]l8|, [HJ. We mention that 
the more complicated helicity amplitudes in ref. [|ll]] were decomposed into supersymmetric and 
non-supersymmetric pieces. The FDH scheme was applied to the supersymmetric pieces in order to 
preserve their supersymmetry; the same functions will thus appear in a supersymmetric decomposi- 
tion of the two-quark-two-gluon and four-quark amplitudes. 

In this paper we also evaluated the two-loop relation between the standard MS coupling and 
the couplings in the four-dimensional helicity and dimensional reduction schemes. It turns out that 
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the couplings in the latter two schemes are identical, at least in the case of a single gauge coupling. 
The two-loop coupling shift was then used to obtain the three-loop gauge theory (5 function in the 
dimensional reduction and four-dimensional helicity schemes, using the previously calculated |5§| , 
59), H MS ^-function. 

It would be interesting to investigate the super symmetry-preserving properties of the FDH 
scheme more generally, including amplitudes with external fermions. Some open issues are whether 
the FDH scheme continues to preserve supersymmetry beyond two loops, and whether there are any 
subtleties analogous to the evanescent couplings that appear in the dimensional reduction scheme. 
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